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ABSTfiACT 


Predicting  officer  attrition  is  a  najor  difficulty  in 
the  Planning,  Progranning,  and  Budgeting  Process.  The 
darine  Corps  has  historically  used  aggregated  estiaates  for 
attrition  rates.  The  purpose  of  this  paper  is  to  demon¬ 
strate  the  application  of  the  Jaaes-Stein  and  other 
shrinkage  type  parameter  estimation  schemes  for  the  purpose 
of  generating  more  stable  manpower  loss  rates.  It  has  been 
illustrated  within  this  project  that  improvement  can  be 
attained  by  application  of  the  James-Stein  and  ninimax 
shrinkage  methods  rather  than  the  more  natural  dazimum 
Likelihood  estimation  process.  The  shrinkage  schemes 
employed  herein  offer  powerful  and  useful  methods  for 
generating  attrition  rates. 
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I.  IHTBODOCTIOS 


A.  GENE6AL 

This  paper  is  the  result  of  a  Headquarters  Harine  Corps 
request  to  explore  the  use  of  James-Stein  and  other 
shrinkage  type  parameter  estimation  schemes  for  the  purpose 
of  generating  manpower  loss  rates.  The  very  large  number  of 
cells  within  the  USHC  officer  force  structure  leads  to  tne 
condition  that  empirical  attrition  rates  are  unstable.  This 
problem  is  compounded  by  the  fact  that  the  cell  probabili¬ 
ties  are  small.  Further  difficulties  are  present  because 
some  of  tne  inventory  cells  are  empty  for  structural  reasons 
while  others  are  empty  by  chance.  Tnerefore,  the  small 
inventory  cells  draw  especial  attention. 

In  order  to  deal  with  this  phenomenonr  variance  stabi¬ 
lizing  and  symmetrizing  transformations  followed  by  versions 
of  the  James-Stein  shrinkage  technique  were  applied  to 
selected  aggregates  of  the  OSHC  officer  data.  The  applica¬ 
tion  of  these  methods  to  the  summary  data  provided  by  Navy 
Personnel  Research  and  Development  Center  (HPRDC)  has  shown 
much  promise.  It  has  been  illustrated  within  this  project 
that  improvement  can  be  attained  from  application  or  these 
methods  over  other  estimating  methods. 

Since  considerable  further  work  is  needed  to  screen  our 
choices,  a  further  purpose  or  this  thesis  is  to  document  a 
number  of  supporting  materials,  including  a  thorougn  discus¬ 
sion  of  the  following: 

1.  Background  of  attrition  and  promotions  for  JSilC 
officers. 

Data  types  and  the  immediately  available  data  rase. 


2. 


3.  Identifying  the  plans  (consumers)  that  utxlize 
attrition  rates. 

4.  Formulas  and  mataematical  risk  calculations  relating 
to  the  James- Stein  and  empirical  Bayes  estimators. 

Cross  validation  procedures  were  utilized  to  illustrate 
the  improvement  of  shrinkage  methods  over  the  use  of 
empirical  rates. 

B.  BiCKGfiOOND 

The  forecasts  made  by  Manpower  planning  models  are 
affected  by  three  general  factors:  existing  inventory, 
projected  losses,  and  projected  gains.  In  order  to  project 
the  inventory  into  various  future  time  periods  it  is  neces¬ 
sary  to  use  a  realistic  system  of  flow  rates.  Currently 
those  rates  are  generated  largely  through  arbitrary  manual 
processes.  The  goal  of  this  paper  is  to  examine  ways  to 
develop  objective  and  computerized  loss  rates. 

There  is  an  inherent  confusion  in  terminology  of  losses 
and  loss  rates  in  that  some  are  leavers  from  the  Marine 
Corps  and  some  are  merely  transitioning  within  the  force 
structure.  Flows  of  Officers  from  particular  cells  are 
Characterized  by  Military  Occupational  Specialty  (MOS) , 
length  of  service  (LCS),  grade,  and  perhaps  some  other  cross 
classifying  characteristics.  When  the  flows  go  from  one 
cell  to  anotner  within  the  Marine  Corps,  then  the  flow  is 
referred  to  as  a  strength  loss  ror  tne  former  ceil.  An 
officer  moving  from  one  cell  tc  anotner  cell  does  not  mean  a 


loss  from  the  Marine  Corps  to  tne  civilian  lacor  market. 
Changes  such  as  promotion  to  a  higher  pay  grade  or  losses  to 
a  cell  due  to  an  officer  increasing  from  one  length  of 
service  to  tne  next  are  other  forms  of  transition  within  tne 
force  structure. 
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losses  from  withiji  the  Narine  Corps  to  the  civilian 
labor  market  are  the  only  losses  treated  in  this  thesis. 
They  can  be  voluntary  or  involuntary.  Voluntary  losses 
occur  when  officers  resign,  retire,  or  are  released  by 
choice.  Involuntary  losses  occur  due  to  discharge,  death, 
disability,  release  from  active  duty,  and  retirement. 
Retirement  can  be  voluntary  or  involuntary. 

Authors  who  have  studied  loss  rates  and  how  they  apply 
to  manpower  planning  are  emphasizing  the  importance  of 
understanding  these  rates  in  order  to  conduct  proper 
forecasting.  [Bef.  1]  [Bef.  2] 

A  difficulty  encountered  in  doing  this  project  was  that 
of  obtaining  sufficient  data  points  with  which  to  conduct 
cohort  or  census  manpower  data  analysis  [Bef.  1].  The 
reason  for  this  difficulty  was  the  state  of  the  Marine  Corps 
automated  data  reporting  system  prior  to  1976.  During  tne 
period  1970  to  197b,  tne  Marine  Corps  was  instituting  its 
first  automated  personnel  reporting  system.  Many  problems 
were  encountered  during  this  period.  Therefore,  data  from 
this  period  is  generally  unreliable.  Data  at  the  Defense 
Manpower  Data  Center  is  based  on  input  from  Headquarters 
Marine  Corps  and  was  found  to  have  similar  problems.  A 
complicating  factor  during  this  period  was  the  Vietnam 
conflict  and  the  attritions  generated  by  the  conflict. 
These  problems  prevented  an  extended  data  base  from  which  to 
draw  from  for  the  purpose  of  this  paper.  The  Marine  Corps 
manpower  reporting  system  was  amended  in  1976  and  by  1977, 
the  system  was  reporting  data  with  over  95  percent  accuracy. 

C.  IBE  MARINE  COBPS  OFFICES  ATTRITION  AND  PBOHOTION 
STBOCTOBE 

An  analysis  of  a  Marine  Corps  officer  as  he  moves 
through  his  career  reveals  key  periods  in  wnich  an  officer 


is  nost  likely  to  leave  the  service.  These  periods  coincide 
vitn  promotion  points  in  the  service.  under  normal  circum¬ 
stances,  unrestricted  2nd  Lieutenants  are  promoted  to  1st 
Lieutenants  after  24  months  in  grade.  1st  Lieutenants  are 
promoted  to  Captains  in  the  Fifth  year  of  commissioned 
service.  Captains  to  Majors  during  the  tenth  year.  Majors  to 
Lieutenant  Colonels  during  the  sixteenth  year  of  service, 
and  Lieutenant  Colonels  to  Colonels  durrng  the  tventy  second 
year  of  service.  An  important  influencing  factor  is  the 
fact  that  retirement  benefits  are  not  achieved  until  the 
twentieth  year  of  total  service. 

1 .  Lieutenants 

Lieutenants  normally  have  a  four  or  five  year 
initial  term  of  service  depending  on  their  source  of  commis¬ 
sion  and  Military  Occupational  Specialty  (NOS)  ,  either  air 
or  ground.  Aviation  Lieutenants  are  some  what  special. 
Only  a  small  numoer  of  aviation  Lieutenants  leave  the 
service.  This  is  due  to  their  longer  initial  term  of 
service,  which  causes  them  to  reach  the  grade  of  Captain 
before  the  level  of  natural  attrition  is  more  likely  to 
occur.  Attrition  of  Aviation  Lieutenants  is  primarily  a 
function  of  accidents,  illness,  and  disciplinary  problems. 

2 .  Captains 

* 

There  is  a  much  greater  period  of  flexibilxty  for 
Captains  given  that  they  are  regular  officers  or  reserves 
witn  an  extension  of  service  agreement.  Captains  are 
promoted  at  five  years  of  commissioned  service  and  are 
either  passed  over  or  promoted  to  Major  by  their  eleventh 
year  cf  service.  Thus  in  a  period  of  over  six  years,  an 
individual  may  leave  the  service  voluntarily  at  any  point. 


D.  AITBIIlOH  OF  AVIATOBS 


The  attrition  rate  of  aviation  officers  was  affected  by 
the  intiation  of  Aviation  Officer  Continuation  Pay  (AOCP)  in 
1981.  AOCP  provides  a  bonus  per  year  to  avxation  officers 
and  this  in  turn  obligates  continued  service.  The  program 
was  applied  to  all  ranks  provided  the  individual  met  certain 
active  duty  flight  status  requirements.  This  action  by  DOD 
has  had  its  desired  effect  upon  the  retention  of  aviation 
officers.  This  accounts,  in  part,  for  the  reduction  in 
attrition  rates  for  Lieutenant  Colonels  in  1983  and  for 
Majors  and  Captains  in  1962  and  1983. 

E.  TERHS 

The  following  terms  will  be  utilized  within  this 
project: 

1.  Aggregation  -  collection  of  historical  officer  data 
over  distinguishing  characteristics.  For  example, 
the  aggregation  of  all  grades  within  a  Military 
Occupational  Specialty  (MOS)  . 

2.  Attrition  -  any  departure  from  the  Marine  Corps  ny  an 
officer. 

3.  Attrition  rate  -  idealized  ratio  of  attritions  to 
inventory. 

4.  End-strength  -  total  number  of  the  Marine  Corps 
officers  at  the  end  of  a  specified  period  of  time. 

5.  Limited  Duty  Officer  (LCO)  -  professionally  qualified 
officers  specifically  designated  for  limited  duty 
within  certain  MOS's. 

6.  Unrestricted  -  all  ofricers  not  designated  as  Limited 
Duty  Officers  (LDO) . 

7.  Regular  -  unrestricted  officers  not  specifically 
assigned  to  a  "Reserve”  status. 
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e  -  unrestricted  officers  with  a  specifically 
ed  "Beserve"  status. 


s\\’ 
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II.  STflgCTgHg  Of  DATA  BASE 


A.  SOHHABI  DATA  FILE 

The  data  used  in  this  project  was  obtained  from  a 
sttmoaiy  data  file  from  the  Commander,  Navy  Personnel 
Research  and  Development  Center  (SPRDC) .  The  summary  data 
file  mas  generated  from  two  source  files:  the  Headquarters 
Master  File  (HHF)  and  the  Quarterly  Statistical  Transaction 
File  (STATS)  .  The  Headquarters  Master  file  (HMF)  was  used 
to  produce  nistorical  officer  inventories  as  of  the  begin¬ 
ning  cf  the  fiscal  year.  Inventories  were  generated  for 
each  fiscal  year  from  1977  through  fiscal  year  1983.  The 
inventories  were  identified  by  distinct  characteristics. 
These  characteristics  were  Military  Occupational  Specialty 
(MOS)  ,  grade  and  length  of  service. 

The  data  in  hand  does  not  distinguish  between  unres¬ 
tricted  officers  and  limited  duty  officers.  For  that 
reason,  discussion  of  the  characteristics  that  separate 
these  two  types  cannot  be  used  advantageously  herein.  Such 
discussion  is  included  however  in  anticipation  of  receiving 
a  more  discriminating  data  base  in  the  future. 

The  Quarterly  Statistical  Transaction  File  was  used  to 
generate  historical  losses.  within  tois  file,  if  the  Type 
Transaction  code  indicated  a  loss,  tnen  the  Effective  Date 
of  Action  field  would  specify  the  year  and  month  of  the 
loss.  Losses  were  classified  into  eight  categories  for  each 
fiscal  year  1977  through  1983.  The  losses  were  furtner 
classified  into  distinct  elements  by  MOS,  grade  and  length 
of  service. 
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B.  SOHBiBX  DATA  FOEflAI 


The  sumaary  data  file  classified  the  Marine  Corps 
Officer  inventory  into  40  Military  Occupation  Specialties, 
10  grade  levels,  31  lengths  of  service,  and  8  loss  catego¬ 
ries.  The  data  format  is  defined  by  Tables  1  and  II.  The 
colunn  containing  the  letters  A  through  £  refer  to  the 
structural  zero  problem  discussed  in  the  next  section.  See 
the  material  surrounding  the  contents  of  Table  V. 

Throughout  the  remainder  of  this  project,  when  reference 
is  made  to  a  particular  Military  Occupational  Specialty 
(MGS)  ,  the  data  code  reference  from  Table  II  will  be  used 
instead  of  the  actual  MOS.  For  example,  this  project  will 
refer  to  the  Utilities  MOS  as  number  06  and  not  11.  It 
should  also  be  understood  that  the  two  digit  NOS  identifier 
listed  in  Table  II  is  strictly  the  Occupational  Field  iden¬ 
tifier  in  the  USMC  MOS  Manual.  This  results  in  MOS,  for  the 
purposes  of  this  paper,  being  an  aggregation  of  what  is 
usually  understood  to  be  an  MOS  cell. 

C.  STBUCTUBAL  ABB  SABBLIIIG  ZEBOES 

There  are  two  kinds  of  zero  inventory  cells  that  must  be 
dealt  with: 

1.  Structural  zeroes;  A  cell  wnose  inventory  is  always 

zero  because  certain  grades  and  length  of  service 
combinations  should  never  appear  in  the  Military 
Occupational  Specialty  (MOS)  .  For  example,  a 

Lieutenant  Colonel  with  4  years  of  service  in  any  MOS 
or  an  infantry  warrant  officer  in  MOS  03  does  not 
exist. 

2.  Sampling  zeroes:  according  to  the  sample  data,  there 
is  no  officer  in  the  particular  grade  and  length  of 
service  combination  for  a  given  fiscal  year  and 
Military  Occupational  Specialty.  This  is  not  a 
permanent  condition,  but  one  tnat  occurs  by  chance. 


TABLE  I 
DATA  FOBHAT 


COLOHN  DATA 

1-2  YEAB  -  (FISCAL)  ( 1  977. .  .  1  983) 

3-4  BOS  GfiOOP  (0...39) 

5-6  GRADE 

7-8  LENGTH  OF  SERVICE  (0...30) 

9-13  BEGINNING  OF  YEAR  INVENTORY 

14-53  LOSS  COUNTS  (TYPES  1...8) 


CODE 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


GRADE 

WARRANT  OFFICER  (W-1) 
CHIEF  WARRANT  OFFICER 
CHIEF  WARRANT  OFFICER 
CHIEF  WARRANT  OFFICER 
SECOND  LIEUTENANT 
FIRST  LIEUTENANT 
CAPTAIN 
MAJOR 

LIEUTENANT  COLONEL 
COLONEL 


(CWO-2) 
(CWO-3 
CWO-4 1 


CODE  TYPE  LOSS 

1  VOLUNTARY  RESIGNATION 

2  VOLUNTARY  REIIREHENT 

3  INVOLUNTARY  -  DEATH 

4  INVOLUNTARY  -  DISCHARGE 

5  INVOLUNTARY  -  DISABILITY 

b  RELEASE  FROM  ACTIVE  DUTY 

7  DISABILITY  RETIREMENT 

8  INVOLUNTARY  RETIREMENT 


For  tne  purpose  of  clarity  in  this  paper,  the  ncnstruc- 
tural  zero  region  will  be  referred  to  as  being  within  the 
Feasible  Region.  This  region  is  defined  by  the  lower  and 
upper  limits  of  the  length  of  service  for  each  grade.  Table 
III  identifies  the  minimum  and  maximum  years  of  commissioned 
service  for  the  Marine  Corps  officer  grades  for  the 


TABLE  II 

HILITAfiX  OCCOPATIOHAL  SPECIALTIES  (HOS) 


DAIA 

CODE 

HOS 

CAT 

MOS  TITLE 

00 

UN 

A 

UNKNOWN 

01 

01 

A 

PERSONNEL  AND  ADMINISTRATION 

02 

02 

A 

INTELLIGENCE 

03 

03 

C 

INFANTRY 

04 

04 

A 

LOGISTICS 

05 

08 

A 

FIELD  ARTILLERY 

06 

11 

D 

UTILITIES 

07 

13 

A 

ENGINEER*  CONSTRUCTION.  AND  EQUIPdENT 
DRAFTING*  SURVEYING,  AND  MAPPING 

08 

14 

0 

OS 

15 

D 

PRINTING  AND  REPRODUCTION 

10 

18 

C 

TANK  AND  AMPHIBIAN  TRACTOR 

1 1 

21 

A 

ORDNANCE 

12 

23 

B 

AMMUNITION  AND  EXPLOSIVE  ORDNANCE 

13 

25 

A 

DISPOSAL 

OPERATIONAL  COMMUNICATIONS 

14 

26 

A 

SIGNALS  INTELLIGENCE/GROUND  ELECTRONIC 

15 

28 

B 

WARFARE 

DATA/COMMUNICATIONS  MAINTENANCE 

16 

30 

A 

SUPPLY  ADMINISTRATION  AND  OPERATIONS 

17 

31 

A 

TRANSPORT  ATI  CN 

18 

33 

A 

FOOD  SERVICE 

19 

34 

A 

AUDITING,  FINANCE,  AND  ACCOUNTING 

MOTOR  TRANSPORT 

20 

35 

A 

21 

40 

A 

DATA  SYSTEMS 

22 

41 

B 

MARINE  CORPS  EXCHANGE 

23 

43 

A 

PUBLIC  AFFAIRS 

24 

44 

A 

LEGAL  SERVICES 

25 

46 

A 

TRAINING  AND  AUDIOVISUAL  SUPPORT 

26 

55 

E 

BAND 

27 

57 

0 

NUCLEAR,  BIOLOGICAL,  AND  CHEMICAL 
MILITARY  POLICE  AND  CORRECTIONS 

28 

58 

A 

29 

59 

B 

ELECTRONICS  MAINTENANCE 

30 

60 

A 

60XX 

31 

61 

A 

AIRCRAFT  MAINTENANCE 

32 

63 

B 

AVIONICS 

33 

65 

B 

AVIATION  ORDNANCE 

34 

68 

B 

WEATHER  SERVICE 

35 

70 

0 

AIRFIELD  SERVICES 

36 

72 

A 

AIR  CONTROL,  AIR  SUPPORT,  AND  ASTI-AIR 

37 

73 

A 

WARFARE 

AIR  TRAFFIC  CONTROL 

38 

75 

C 

PILOTS  AND  NAVAL  FLIGHT  OFFICERS 

39 

99 

E 

IDENTIFYING  HOS  AND  REPORTING  HOS 

TABLE  Ill 


flIHIBOH  AID  HAXIHOE  YEAES-OF-COEHISSIOMED 
SEE VICE  FOE  OFFICEE  GEADES 


Grade 

Unrestricted 

HIN  yes  MAX  yes 

LOO 

MiE  yes 

MAX  yes 

HOI 

1 

4 

CW02 

4 

9 

— 

— 

CiJ03 

8 

13 

— 

— 

CU04 

12 

24 

— 

— 

21T 

1 

4 

1LI 

1 

8 

1 

8 

CAPT 

4 

13 

2 

12 

HAJ 

9 

20 

5 

18 

IICOL 

14 

23 

8 

23 

COL 

21 

30 

— 

- 

In  order  to  get  the  o£ficer*s  length  of  service  one  must 
add  the  officer's  prior  enlistment  period  to  the  officer's 
years  of  commissioned  service.  Since  prospective  warrant 
officers  are  reguired  to  have  a  minimum  of  five  years  of 
service  prior  to  consideration  by  the  officer  selection 
board,  the  logical  lover  limit  for  length  of  service  for 
warrant  officers  is  five  years.  The  LOS  then  would  include 
the  total  of  this  service  and  the  respective  years  of 
commissioned  service.  Additionally,  warrant  officers  are 
reguired  to  have  10  years  of  service  prior  to  consideration 
to  become  Limited  Duty  Officers.  This  affects  the  lower 
limit  of  the  length  of  service  for  Limited  Duty  Officers. 
Tae  logical  lower  limit  for  length  of  service  for  LDO's 
wouxd  be  10  years. 

The  summary  data  provided  for  this  project  did  not 
distinguish  between  unrestricted  and  Limited  Duty  Officers. 
Additionally,  the  length  of  service  was  the  characteristic 
provided  instead  of  the  years  of  commissioned  service  for 
each  of  the  officer  grades.  Table  IV  identifies  the  lower 


and  upper  limits  of  the  Feasible  Region  of  the  length  of 
service  for  Narine  corps  officer  grades  utilized  in  this 
project. 
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TABLE  IT 

NIIIHON 

AND  HAXIHOH  LENGTH  OF 

SEBTICE 

FOB  OFFICEfi  GBABES 

1 

Grade 

HIM  LOS 

NAX  LOS 

HOI 

5 

20 

CH02 

8 

30 

CH03 

12 

30 

CUC4 

16 

30 

2LT 

1 

20 

UT 

1 

20 

CAPT 

4 

30 

NAJ 

9 

30 

LTCOL 

14 

30 

COL 

1 

21 

30 

Mote  the  maximum  limit  in  the  majority  of  grades  can 
feasibly  be  30  years.  This  was  substantiated  by  the  actual 
data  wherein  the  maximum  length  of  service  for  grades  CW03 
through  Colonel  was  30  years.  The  maximum  length  of  service 
for  hCI's,  2LT's  and  iLT’s  was  limited  to  20  years  because 
of  the  very  small  number  of  actual  officers  with  21  to  30 
years  of  service  within  these  grades. 

One  method  of  considering  the  size  of  the  feasible 
region  is  to  create  a  multidimensional  array  using  all 
combinations  of  different  characteristics  within  the  data. 
If  ail  characteristics  were  considered  for  all  seven  fiscal 
years,  the  size  of  the  multidimensional  matrix  wculd  be 
7x40x31x10x8.  This  would  represent  7  years,  40  MOS's,  31 
LOS's,  10  grades  and  8  loss  types.  This  would  create 
694,400  different  cells  that  could  be  dedicated  to  loss 


counts.  However,  to  bring  the  dimension  to  a  more  manage¬ 
able  level,  consider  one  year  o£  data  for  each  loss  type 
having  a  cell  count  of  40x31x10.  Therefore,  for  each  loss 
type  within  a  fiscal  year  there  is  a  total  of  12,400  cells. 
Hany  cf  these  cells  are  infeasible  and  are  designated  those 
with  structural  zeroes. 

The  region  of  structural  zeroes  is  the  area  outside  of 
the  minimum  and  maximum  length  of  service  for  the  officer 
grades  as  listed  in  Table  IV.  But  it  is  further  articulated 
because  of  additional  constraints  imposed  by  the  various  MOS 
classifications.  This  is  largely  because  of  the  suitability 
of  certain  MOS*s  being  assigned  to  unrestricted,  warrant  and 
Limited  Duty  Officers.  The  particular  dOS*s  are  identified 
for  illustration  purposes  by  letter  in  Tanle  ll  according  to 
the  general  categorization  listed  in  Table  V.  This  general 
aos  categorization  specifies  by  category  the  inclusive 
grades  within  certain  MOS*s,  the  number  of  MOS's  with  the 
specified  grade  structure,  the  number  of  structural  zeroes 
per  hOS,  and  the  total  structural  zeroes  within  each 
category.  The  total  number  of 'structural  zeroes  is  6149. 
Therefore,  the  number  of  cells  within  the  feasible  region  is 
625  1  for  each  loss  type  for  a  given  fiscal  year. 

D.  C£liTfiAL  DATA 

The  following  is  an  observation  made  by  the  author 
regarding  the  summary  data  utilized  in  this  project.  This 
data,  which  was  provided  by  Navy  Personnel  Research  and 
Development  Center  (NPROC),  was  classified  as  central  data 
according  to  Bartholomew  [Ref.  1:  p.  25].  Bartholomew 

defines  central  wastage  (attrition)  rates  as  follows: 

1.  The  central  wastage  (attrition)  rate  is  the  number  of 
leavers  during  the  period  who  were  in  this  class  when 
they  left  divided  by  the  average  number  in  this  class 
during  the  period. 


TABLE  ? 

SIBUCTOBAL  ZEBCES  CATEGOBIES 


Grades 

Number 

Stru. 

Zeroes 

Total 

Zeroes 

Category 

Within  MOS 

of  MOS 

per  MOS 

per  Cat. 

A 

W01... IICOL 

23 

129 

2967 

B 

W01...  CW04,LD0 
2LI...  LICOL 

8 

159 

1272 

C 

3 

202 

606 

D 

U01  ...  CU04 

5 

237 

1185 

E 

W01...CCL 

1 

119 

119 

TOTAL 

40 

6149 

j 

A  problem  arose  on  several  occasions  when  the  data  was 
disaggregated  to  a  level  where  the  inventory  was  very  small. 
For  example,  when  examining  the  inventory  in  a  particular 
fiscal  year,  the  inventory  was  zero  ror  a  Length  of  Service 
(A-OS)  and  Military  Occupational  Specialty  (MOS)  combination. 
Examining  the  inventory  in  the  next  fiscal  year  for  the  same 
LOS  and  MOS  combination  also  was  zero.  However,  the  problem 
arises  when  the  numner  of  leavers  is  equal  to  or  greater 
than  cne.  The  average  of  the  inventories  for  the  two  fiscal 
years  is  zero.  Using  this  result  in  tne  estimation  of  the 
attrition  rate  would  be  ambiguous.  The  zero  in  the  denomi¬ 
nator  of  the  ratio  of  leavers  to  average  inventory  resalts 
in  an  undefined  expression. 

For  the  purpose  of  removing  this  ambiguity  in  the  data, 
the  following  policy  was  adopted  to  define  the  central 
inventory  number  for  the  officer  force  at  disaggregated 
levels:  For  any  cell  or  collection  of  cells 

1.  let,  t  =  1...6,  refer  to  the  years  1 977.  . .  1  i?a2. 

2.  Let,  y  (t)  =  Number  of  losses  in  year  t. 

3.  net,  INV(t)  =  Inventory  in  the  beginning  of  year  t. 
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.  Let,  N  (t)  =  Haximam  of  y  (t)  and  the  average  inventory 
in  year  t,  using  the  beginning  inventory  in  year  t 
and  t+l  and  computing  their  average: 

(INV(t)  +INV(t  +  l))/2. 

5.  Let,  m  =  Sum  of  ¥(t)  divided  by  sum  of  N(t),  (both 
sums  over  t)  represents  the  empirical  central 
attrition  rates  for  the  particular  cell  or  aggregate. 

The  data  bank  that  will  ultimately  be  utilized  in  the 
parent  project  at  NPfiCC,  of  which  this  project  is  a  subset, 
will  contain  the  appropriate  inventory  and  attrition  data  so 
that  the  above  modifications  will  not  be  required.  It  is 
currently  planned,  hy  NPRDC,  to  produce  the  historical 
officer  inventory  as  of  the  beginning  of  a  month  and  tne 
losses  will  be  identified  by  the  month  in  which  the  officer 
leaves  the  Marine  Corps.  This  will  provide  the  elements  for 
a  more  accurate  estimation  of  tne  attrition  rate  on  the 


III.  COBHEMT  SYSTEM 


A.  OYEBVIEH 

In  order  to  gain  an  appreciation  for  the  complexity  of 
the  Manpower  Planning  process,  the  author  during  his  experi¬ 
ence  tour  at  Headquarters  Marine  Corps  had  the  opportunity 
to  contact  key  individuals  responsible  for  the  execution  of 
tne  models  currently  utilized  within  the  Manpower  depart¬ 
ment.  Ihe  methodology  used  to  identify  the  current  system 
included  interviews  with  the  action  officers  of  the  Officers 
Plans  Section  (Code  MPP-30)  ,  briefings  and  interviews  .with 
representatives  from  Allocations  Section  of  Manpower  Control 
Branch  (Code  MPC-20)  ,  Officer  Assignment  Branch  (Code  MMOA)  , 
Manpower  Plans,  Programs  and  Budget  Section  (Code  MPP-40) , 
Manpower  Management  Information  Systems  Branch  (Code  MPI) , 
Assistant  Deputy  Chief  of  Staff  for  Manpower  (Code  M) ,  and 
examination  of  existing  manuals,  Beguirements  Statement  for 
Officer  Planning  System,  and  Statement  of  Work  for  the 
Develcpment  of  the  Officer  Inventory  and  Manpower  Flow  Data 
Base  System. 

Ibe  goal  of  the  officer  planning  process  is  to  shape  the 
officer  force  structure  over  the  7  planning  years  (current, 
budget,  and  five  Program  Objective  Memorandum  (POM)  years)  . 
This  fence  structure  is  affected  by  three  influencing 
factors:  existing  inventory,  projectea  losses,  and  expected 
gains.  The  Officer  Plans  Section  (Code  MPP-30)  supports  the 
Manpower  Department  at  Headquarters  Marine  Corps  by 
preparing  plans  that  address  each  of  these  factors. 

The  plans  which  input  the  target  inventory  calculations 
are  the  Promotion  Plan,  the  Lateral  Move  and  Directed 
Lateral  Move  Plans,  the  Inventory  Projection,  and  the 
Manpower  Plan. 


Losses  la  the  officer  force  are  the  result  of  retire¬ 
ments,  resignations,  releases,  discharges  and  other  factors 
(e. g.,  death).  Among  these  factors,  releases  and  mandatory 
(statutory)  retirements  addressed  in  the  Oficer  Inventory 
Control  Plan,  can  he  estimatea  with  a  high  degree  of 
certainty.  The  remaining  loss  factors  are  estimations  based 
on  forecasted  attrition  rates  utilizing  historical  data. 
All  of  the  factors  contributing  to  losses  in  the  officer 
inventory  are  addressed  in  the  Inventory  Projection  Loss 
Plan. 

Gains  in  officer  force  are  tne  result  of  accessions  and 
augmentation.  Accessions  are  identified  in  the  Accession 
Plan  and  augmentations  are  identified  in  the  Augmentation 
Plan.  The  number  of  field  grade  officers  are  directly 
influenced  by  statutory  limitations  specified  in  Title  10, 
(J.S.  Code.  These  limitations  affect  both  of  these  plans. 
In  the  planning  process,  in  order  to  plan  for  an  adequate 
number  of  field  grade  officers  in  the  future,  one  must 
currently  access  a  sufficient  number  of  second  lieutenants. 

The  force  structure  is  also  affected  by  the  mixture  of 
skills  of  the  officer  inventory.  The  plans  used  to  support 
the  required  mix  and  levels  of  skills  are  the  TSS/dOS 
Distribution  Plan,  the  Entry  Level  Training  Plan,  tne 
lateral  Mover  Plan,  the  Continuation  Plan,  and  the  Promotion 
Plan . 

B.  SCUfiCES  OF  IlFOBHAIION 

Input  from  the  following  external  sources  provide  infor¬ 
mation  required  by  the  Officer  Plans  Section  to  create  the 
plans  identified  above. 

1.  Statutory  limitations  specified  in  Title  10,  J.S. 

Code. 
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2.  Desired  iaventcry  force  structure  as  enumerated  in 
the  Grade  Adjusted  Recapitulation  (GAB). 

3.  FeedhacJc  data  from  agencies  which  execute  the  plan 
(for  example,  actual  monthly  promotions) . 

4.  Historical  data  from  JUMPS/MMS  on  officer  service 
characteristics  and  any  changes  in  the  officer 
inventory  over  time. 

5.  Current  data  on  the  existing  officer  inventory 
utilizing  the  Officer  Slate  file  and  the  Lineal  List 
file. 

C.  flON  OF  IHFOfiBATIGIl 

Ihe  planning  process  consists  of  three  major  elements 
comprised  of  several  processes.  These  elements  are: 

1.  Preparation  for  planning. 

2.  Plan  development. 

3.  Force  structure  analysis. 

The  preparation  process  includes  development  of  the 
targeted  force  structure  to  include  distribution  of  officers 
by  skill,  sex,  uncestricted/restricted  categories,  and 
grades  based  on  historical  flow  rates  and  USdC  policy  objec¬ 
tives.  The  requirements  definition  process  establisnes  the 
authorized  number  of  officers  in  eacn  military  Occupational 
Specialty  (HOS)  and  grade  for  the  five  planning  years  and 
determines  the  target  force.  The  programmed  requirement  is 
specified  in  the  Numerically  Adjusted  Recapitulation  (NAR)  . 

The  planning  process  also  includes  toe  determination  of 
officer  attrition  rates.  Currently  the  attrition  rate 
determination  method  is  accomplished  in  conjunction  with  the 
Officer  Promotion  Planning  Process.  Inis  model  divides  the 
attrition  rates  into  two  categories:  statutory  and  non- 

statutory.  The  statutory  attrition  rate  results  from  'JS;iC 
and  EOD  policies  and  is  applied  directly  to  the  officer 
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iavestory.  The  non- statutory  rate  is  derived  from  a  time 
series  analysis  of  officer  inventory.  This  rate  is  based  on 
historical  data  and  assumes  ccntinuation  of  present  trends. 
Currently  there  is  no  attempt  to  provide  predictive  rates 
based  on  econometric  techniques.  However,  such  an  approach 
is  being  researched  by  the  Navy  Personnel  Research  and 
Development  Center  for  use  by  the  Manpower  Department. 

The  development  process  includes  projection  of  the 
officer  inventory  and  generation  of  various  plans.  The  goal 
of  the  inventory  projection  process  is  to  provide  a  break¬ 
down  by  month  of  the  number  of  officers  in  each  grade. 
Military  Occupational  Specialty  (MOS) ,  and  other  officer 
attributes  (for  example,  restricted,  years  of  commissioned 
service,  sex,  etc.).  This  breakdown  is  for  each  of  the 
planning  years.  Currently  the  inventory  projection  is 
determined  in  one  of  several  ways,  depending  on  the  plan 
being  developed.  The  Promotion  Planning  Process  projection 
is  utilized  if  the  level  of  aggregation  is  adequate  for  the 
plan.  Otherwise,  manual  calculations  are  used  to  aggregate 
from  the  promotion  planning  process  inventory  projection  or 
to  generate  a  new  projection. 

The  plan  generation  process  is  complicated  by  the  inter¬ 
dependence  of  the  plans  and  the  timing  of  their  production. 
Table  VI  identifies  the  plan  generation  process  including 
each  plan,  the  frequency  of  the  generation  and  update  of 
each  plan,  the  users,  and  the  sources  of  data. 

The  force  structure  analysis  process  involves  verifica¬ 
tion  of  the  generated  plans  and  evaluation  of  their  quality. 
Comparison  is  made  of  the  results  of  the  planning  process 
with  the  USMC  objective  force  specified  by  the  Grade 
Adjusted  Recapitulation  (GAB)  and  the  promotion  flow  points 
specified  by  Defense  Officer  Personnel  Management  Act 
(DOPMA)  . 


TABLE  ?I 

OFFICE B  PLAN  GENEBAIIOH  SOMMABI 


SOOBCES 


PLAN 

GENERATED 

UPDATED 

USED  BY 

OF  DATA 

Promotion  plan 

Annual 

Monthly 

Prom.  Brd 
MPP-30 

RES 

MMPB 

1j2,3,«, 

Augmentation 

Plan 

Annual 

Monthly 

MPP-30 

MMOA 

OBB 

1^2, 3, 4, 

Inventory 

Projection 

Annual 

Annual 

MPP-30 

1 »2, 3, 4, 

Inventory  Proj. 
Loss  Plan 

Annual 

Annual 

MPP-30 

1,2,7 

Continuation 

Plan 

Annual 

Annual 

IlPP-30 

1,2, 4, 5 

Accession  Plan 

Annual 

Monthly 

MPP-30 

MMBD 

Entry  Level 
Training  Plan 

Annual 

Monthly 

MMOA 

TBS 

T 

1,2,4 

TBS/aCS 

Distribution 

Plan 

Annual 

Monthly 

MMOA 

IBS 

m 

X 

1,2,4 

Lateral  Move 

Plan 

Annual 

Monthly 

MPP-30 

MMOA 

1,2, 4, 5 

Manpower  Plan 

Annual 

Monthly 

MPP-30 

MPP-40 

1  , 2 , 3 , 4, 

1),6,8,9 

Offic.er  control 
Plan 

Annual 

Monthly 

MPP-30 

1,2,4, 5 

Source  of  Data  Codes 


1 .  HdF 

2.  Attrition  Bates 

3.  DOPMA 

4.  GAB 

5.  Budget 

b.  Candidate  List 


7.  JUflPS/MtIS 

8.  MO  Accessions 

9.  Retirement/Separation 

10.  Lineal  List 

11.  Oificer  Slate  File 


D.  FIAMS 


Fiofflotion  plan  is  generated  by  the  use  of  the  Promotion 
Planning  Process  model.  The  model's  output  is  a  statistical 
projection  of  the  officer  inventory  over  a  period  of  up  to 
10  years.  It  also  provides  several  built  in  promotion 
strategies. 

Augmentation  Plan  is  generated  to  provide  the  Officer 
Retention  Board  (OBB)  with  information  on  the  vacancies 
available  in  each  fiscal  year  of  the  planning  process 
regarding  officers  in  a  category  (i.e. ,  Aviation,  Ground, 
and  Naval  Flight  Officer)  with  a  certain  number  of  years  of 
commissioned  service. 

Inventory  Projection  Plan  is  generated  for  a  five  year 
projection  and  is  specified  for  each  year  by  grade,  years  of 
commissioned  service,  and  category  (unrestricted  or  Limited 
Duty  Officer  (LDO)). 

Inventory  Projection  Loss  Flan  contains  data  on  expected 
losses  from  the  inventory  projection  process.  The  plan  is 
utilized  to  indicate  which  year  of  commissioned  service  the 
losses  are  most  likely  to  occur. 

officer  Inventory  Control  Plan  provides  information  on 
the  number  of  reserve  officers  whose  contracts  expires 
during  the  fiscal  year.  It  rs  influenced  by  the 
Augmentation  Plan  and  provides  guidance  on  the  number  of  the 
officers  who  may  be  retained. 

Continuation  Plan  specifies  the  number  of  officers  tnat 
may  be  granted  an  extension  of  active  duty  depending  upon 
their  Military  Occupational  Specialty  (M03) . 

Accession  Plan  identifies  the  desired  number  of  incoming 
potential  officers  and  receives  as  input  the  number  of 
vacancies  for  each  year  of  the  planning  process.  Since 
accessions  into  the  USWC  are  by  officer  candidates  who  have 
usually  signed  contracts  several  years  prior  to  graduation. 


the  pcol  of  accessions  for  the  first  three  years  are  gener¬ 
ally  well  known.  The  goals  for  the  remaining  planning  years 
are  established  on  projected  officer  losses  and  on  budgeted 
end  strengths. 

Entry  Level  Training  Plan  and  the  TBS/HOS  Distribution 
Plan  ensure  that  the  entry  level  officers  have  the  mix  of 
Military  Occupational  Specialties  that  is  needed  to  fully 
support  the  force  structure. 

Lateral  Move  plan  provides  information  on  the  shortage 
or  overage  in  each  Military  Occupational  Specialty  by 
year-group. 

Manpower  Plan  is  a  summary  of  gains  and  losses  in  each 
grade . 


IT.  HAHIME  CQgPS  OFFICER  ATTBITIOH  MODEL 


A.  DFFIMIIIOM 

A  deiinitioa  for  a  aaapower  model  offered  Ly  Bartholomew 

is: 

A  manpower  model  is  a  mathematical  description  of  how 
cnaoge  taAes  place  in  the  system.  First  of  all  this 
requires  the  specification  of  any  constraints  under 
wiiich  the  system  operates  ....  Secondly,  a  model 
must  specify  the  mechanism  which  generates  flows.  Some 
flows,  such  as  premotion  or  demotion,  are  under  the 
direct  control  of  management  ....  Other  flows,  such 
as  voluntary  wastage,  are  not  under  direct  management 
control  and  assumptions  a^cut  their  future  levels  is 
likely  to  he  based  on  a  blend  of  historical  data  and 
management  judgement.  £Eef.  1:  p.  7] 


B.  ASSaMPTIOIS 

Assumptions  regarding  flows  of  personnel  within  a  system 
can  be  classified  in  various  ways.  If  a  manpower  policy 
dictates  that  a  specific  percentage  of  individuals  in  a 
particular  category  would  leave  in  a  given  time  period  then 
the  assumption  about  the  flow  would  be  deterministic.  Given 
a  level  of  inventory  in  the  category  under  the  deterministic 
assumption,  then  there  is  no  uncertainty  about  how  many 
individuals  will  actually  leave.  A  deterministic  model 
contains  no  random  variables  and  there  is  a  unique  set  of 
model  outpur  data  for  a  given  set  of  inputs.  [  Bef .  3:  p.  3] 

However,  if  each  individual  in  a  category  had  a  certain 
probability  of  leaving  in  a  given  time  period,  then  the 
assumption  about  the  flow  would  be  stochastic.  Given  a 
level  of  inventory  in  tne  category  under  the  stochastic 
assumption,  then  a  prediction  of  the  number  of  individuals 


leavxiig  would  not  be  precise.  A  stochastic  model  contains 
one  OE  more  random  variables.  The  output  of  a  stochastic 
model  is  random  and  thus  only  estimates  the  true  character¬ 
istics  of  the  system  [Bef.  3:  p.  3].  The  output  of  a 
stochastic  model  is  a  probability  distribution.  If  individ¬ 
uals  within  a  system  behave  independently  regarding  their 
decision  of  leaving,  the  actual  number  of  leavers  has  a 
bxnomial  distribution  whose  average  is  a  certain  percentage 
of  the  inventory  within  a  category  [Ref.  1:  p.  7]. 
Voluntary  attrition  from  the  Marine  Corps  is  a  stochastic 
flow  since  it  results  from  a  number  of  more  or  less  indepen¬ 
dent  individual  decisions. 

C.  A 1TB HIGH  BATES 

A  fundamental  role  in  manpower  analyses  is  played  by 
attrition  rates.  The  time  period  can  be  any  duration.  The 
most  ccmmonly  used  periods  are  months  and  years.  The  use  of 
the  central  attrition  rate  when  the  categories  are  defined 
with  respect  to  length  of  service  can  be  considered  as  esti¬ 
mating  the  constant  rate  of  leaving  or  separation  during  a 
year.  Using  this  method,  the  standard  error  of  the  central 
attrition  cate  can  be  estimated  if  some  assumptions  are  made 
about  the  attrition  process: 

1.  The  number  in  the  category  throughout  the  year  was 
constant  and  equal  to  the  average  number  in  txiis 
class  during  the  year. 

2.  Let,  m  =  the  central  attrition  rate. 

3.  Let,  S*  =  Average  number  in  a  category  during  the 
year. 

4.  Let,  L*  =  Number  of  leavers  during  the  year  who  were 
in  this  category  when  they  left. 

5.  Then,  m  =  L*  /  S*. 


6.  Each  individual  in  the  category  is  subject  to  the 
probability  (mdz)  o£  leaving  in  each  interval 
(z,z-*-dz).  tBef.  1:  p.  25] 

Under  these  assuaptions  the  number  of  losses  would  be 
distributed  Poisson  (3*01).  The  estimated  standard  deviation 
or  the  central  attrition  rate  would  be:  [Ref.  1:  p.  25] 

sd  (m)  =  m  /  sgrt  (L*) 

D.  HARIIE  CORPS  OFFICER  CENTRAL  ATTRITION  RATES 

Seme  eopirical  attrition  rates  are  provided  for  the 
purpose  of  familiarization.  First  are  macro  (large  aggre¬ 
gate)  examples  and  proceed  progressively  to  more  refined 
categorized  examples.  An  example  of  the  estimated  central 
attrition  rates  and  their  estimated  standard  deviations  are 
provided  in  Table  VII.  The  average  inventory  over  the  seven 
year  period  is  also  provided  in  Table  VII  in  order  to 
illustrate  the  number  of  officers  within  each  MOS. 

The  rates  in  Table  vii  are  an  aggregation  of  all  grades 
within  a  Military  Occupational  Specialty  (HOS)  and  all  loss 
types.  The  time  period  includes  all  seven  years  of  inven¬ 
tory  counts  and  the  losses  for  fiscal  years  1977  to  1982. 
The  central  attrition  rate  was  calculated  as  follows: 

1.  Let,  t  =  1...6,  covering  years  1977. ..1982. 

2.  Let,  I  (t)  =  Number  of  losses  in  year  t. 

3.  Let,  lNV{t)  =  Inventory  in  the  beginning  of  year  t. 

4.  Let,  N (t)  =  Average  inventory  in  year  t,  using  the 
neginning  inventory  in  year  t  and  t+1  and  computing 
their  average:  (INV  (t)  ♦INV  (t+ 1)  ) /2. 

5.  Let,  m  =  Sum  of  all  Y(t)  divided  by  sum  of  all  N(t). 

The  central  attrition  for  each  Military  Occupational 

Specialty  is  illustrated  for  the  length  of  service  period 
from  0  to  30  in  Appendix  A,  in  figures  A. 1  toA.10.  Tnese 
figures  illustrate  tne  aggregate  central  attrition  rate  for 


TABLE  7II 

CEHTBAL  ATTBIIlOH  RATES  AND  STANDARD  DEVIATIONS 


NOS 

AVERAGE 

INVENICRY 

BATES 

STANDARD 

DEVIATIONS 

01 

531 

.  121 

.006 

02 

288 

.119 

.00  8 

03 

3059 

.091 

.00  2 

04 

183 

.104 

.010 

05 

1322 

.  101 

.004 

06 

36 

.080 

.019 

07 

657 

.114 

.005 

08 

c 

.034 

.034 

09 

6 

.027 

.027 

10 

512 

.  105 

.006 

1 1 

102 

.090 

.012 

12 

76 

.  100 

.015 

13 

782 

.  138 

.005 

14 

1  28 

.099 

.01  1 

15 

94 

.098 

.013 

16 

11  48 

.  139 

.004 

17 

26 

.  156 

.032 

18 

40 

.  120 

.022 

19 

309 

.  109 

.008 

20 

364 

.115 

.008 

21 

251 

.  143 

.010 

22 

22 

.156 

.035 

23 

59 

.  150 

.020 

24 

459 

.  193 

.008 

25 

24 

.  124 

.02  9 

26 

17 

.081 

.029 

27 

1 1 

.047 

.033 

28 

156 

.093 

.010 

29 

119 

.  120 

.013 

30 

262 

.100 

.008 

31 

0 

.000 

.000 

32 

65 

.042 

.017 

33 

69 

.067 

.01  1 

34 

14 

.108 

.036 

35 

32 

.088 

.021 

36 

526 

.  125 

.006 

37 

108 

.091 

.012 

38 

5644 

.081 

.002 

39 

1296 

.  100 

.004 

TOTAL 

18780 

.  102 

.00  1 

all  eight  types  of  losses  over  the  sxx  years  froB  1977  to 
1982.  These  figures  demonstrate  the  breakpoints  and  discon¬ 
tinuities  within  each  NOS  over  the  distribution  of  length  of 
service.  One  can  very  easily  note  the  increased  rates  wnich 


occur  at  certain  years  of  length  of  service.  The  increased 
rates  frequently  occur  at  the  four,  20,  and  30  years  of 
length  of  service.  The  central  attrition  rate  for  each  iHOS 
is  illustrated  for  the  grades  from  Warrant  Officer  to 
Colonel  (labelled  0. ..9  on  the  abscissa)  in  Appendix  A,  in 
Figures  A. 11  to  A. 20.  These  figures  also  illustrate  the 
aggregate  central  attrition  rate  for  all  eight  types  of 
losses  over  the  same  six  years. 

The  central  attrition  rate  for  each  fiscal  year  aggre¬ 
gated  ever  all  categories  is  provided  in  Table  VIII  and 
illustrated  in  Figure  4.1.  The  trend  is  stable  for  the 
first  three  years  and  a  decreasing  trend  is  evidenced  in  the 
last  three  years. 


1 - 

OSMC  OFFICEB 

TABLE  VIII 

LOSS  BATES  PEB 

FISCAL 

—  1 

lEAB 

1977 

1978 

1978 

1980 

1981 

1982 

.1133 

.  1137 

.1136 

.  1012 

.0928 

.0781 

J 

Table  IX  is  provided  in  order  to  gain  an  appreciation  of 
the  number  of  losses  within  each  loss  type.  This  table 
includes  the  number  of  losses  for  each  fiscal  year  1977 
through  1983  aggregated  over  all  grades,  LOS,  and  MOS.  This 
table  also  identifies  the  loss  types.  One  can  note  the 
large  numner  of  losses  throughout  the  seven  years  in  types 
two,  four,  and  six.  Additionally,  atter  reaching  a  total 

the  general  trend  is  a  decreasing  number  of 


higii  in  1978, 


USMC  OFFICER  LOSS  RATES  PER  FISCAL  YEAR 


FISCAL  YEARS 


Figure  4.1  OSflC  Officer  Less  Bates  Per  Fiscal  Year. 


The  central  attrition  rates  tor  this  level  or  aggregation 
lor  fiscal  years  1977  tnrough  1982  arc  iiiustiatea  in  Taijle 


The  central  attrition  rates  for  each  Military 
Occapaticnal  Specialty  (NOS)  ana  for  each  tiscal  year  l':»77 
through  1982  is  proviUeu  in  Tat-le  II.  These  rates  are 
aggregated  over  all  less  types.  These  rates  are  graphically 
illustrated  in  Appendix  A,  in  Figures  A. 21  to  A. 2b.  one  can 
observe  the  change  of  rates  over  time  tor  each  Nt^S.  It  can 


T&BLE  IX 

LOSSES  PEB  FISCAL  YEAR  BZ  LOSS  TYPE 


LOSS 

TYPE 

1977 

1978 

1979 

1  980 

1981 

1982 

1983 

TOTAL 

11 

118 

90 

109 

40 

24 

37 

35 

453 

L2 

675 

701 

681 

594 

531 

524 

362 

4068 

L3 

32 

42 

38 

39 

34 

22 

27 

234 

LU 

364 

515 

53  6 

572 

532 

431 

317 

3267 

15 

1 

1 

0 

0 

1 

0 

0 

3 

L6 

894 

772 

750 

595 

522 

406 

436 

4375 

17 

15 

27 

17 

2 

3 

2 

10 

76 

18 

52 

50 

39 

35 

43 

30 

4  1 

290 

TOTAL 

215  1 

21S8 

2170 

1877 

1690 

1452 

1228 

12766 

LOSS  TYPES 

Voluntary  Hesignation 
Voluntary  Retirement 
Involuntary  -  Death 
Involuntary  Disccarge 
Involuntary  Disability 
Release  trom  Active  Duty 
Disability  Retirement 
Involuntary  Retirement 


be  noted  with  interest  that  the  MOS's  with  larger  invento¬ 
ries  (refer  to  Table  VII  for  the  average  inventory  witnin 
each  fiCS)  do  not  change  as  much  over  time  as  the  MCS's  with 
smaller  inventories.  For  example,  the  infantry  KOS  OJ  has 
tJie  second  largest  inventory  count.  The  range  of  the  rates 
over  time  fluctuate  from  a  low  of  0.070  to  a  high  of  0.106 
giving  a  range  of  0.036.  The  logistic  dOS  04,  on  the  otner 
nani,  has  a  much  smaller  inventory  count  and  the  rates  over 
time  fluctuate  from  a  low  of  0.049  to  a  high  of  0.22,  giving 
a  range  of  0.171.  In  ether  words,  tne  total  inventory  count 
in  a  category  nas  a  definite  effect  on  the  change  of  rates 
over  time. 


TABL£  X 

CEMTBAL  ATTfilTIOH  BATE  PEB  FISCAL  IBAB  BY  LOSS  TYPE 


LOSS 

TYPE 

1977 

1978 

1979 

1980 

1981 

1982 

LI 

.006 

.005 

.  006 

.002 

.001 

.002 

L2 

.03  6 

.036 

.  035 

.031 

.028 

.028 

L3 

.00  2 

.002 

.  002 

.002 

.002 

.001 

L4 

.019 

.027 

.  028 

.030 

.028 

.023 

L5 

.00  0 

.000 

.  000 

.000 

.000 

.000 

L6 

.047 

.040 

.  039 

.031 

.027 

.022 

L7 

.00  1 

.001 

.  001 

.900 

.000 

.000 

L8 

.00  3 

.003 

.  002 

.002 

.002 

.002 

The  change  in  the  rates  are  also  arfected  ty  the  level 
of  disaggregation.  The  central  attrition  rates  for  the 
major  military  Occupational  Specialty  groups  (Aviation, 
Comnat  Support,  Ground  Combat  and  total)  are  listed  in  Table 
XII.  These  major  groups  are  illustrated  in  Appendix  A,  in 
Figures  A.2S,  A. 30,  A. 31,  and  A. 32.  The  rates  demonstrate 
the  aggregation  over  all  types  of  losses.  For  brevity  tne 
illustrations  were  limited  to  the  grades  ox  Second 
Lieutenant  to  Lieutenant  Colonel. 

A  further  disaggregation  is  listed  in  Table  XIII  and 
illustrated  in  Figure  4.2.  These  rates  demcnstrate  the 
aggregation  within  the  Combat  Support  group  over  all  lengths 
of  service  (LOS)  for  the  voluntary  retirements  (type  2  loss) 
for  the  grades  of  First  Lieutenant,  Captain,  Major  and 
lieutenant  Colonel.  A  common  oiservation  within  this  loss 
type  is  the  increased  rate  as  tne  grade  increases.  This 
pheaciienon  also  corresponds  to  the  increased  rate  of 
attrition  in  the  larger  LOS  categories. 

Another  level  of  disaggregation  of  the  Combat  Support 
group  is  listed  in  Table  XIV  and  illustrated  in  Figure  a. 3. 
These  rates  demonstrate  the  aggregation  within  tne  Motor 


TABLE  Z1 

CENTBAL  ATTBITIOli  BATES  FOR  FISCAL  YEARS  BX  HOS 


cs 

1977 

1978 

1979 

1980 

1981 

1982 

1 

.  132 

.  137 

.  134 

.  132 

.105 

.092 

2 

.  136 

.  151 

.111 

.  112 

.120 

.093 

3 

.  106 

.092 

.098 

.  088 

.086 

.070 

4 

.049 

.  220 

.079 

.097 

.128 

.085 

5 

.110 

.  103 

.  120 

.118 

.079 

.073 

6 

.057 

.  121 

.  131 

.029 

.133 

.024 

7 

.  128 

.  119 

.113 

.  109 

.124 

.087 

8 

.  286 

.000 

.00  0 

.000 

.000 

.000 

S 

.000 

.000 

.200 

.000 

.000 

.000 

10 

.  135 

.097 

.  150 

.094 

.091 

.058 

11 

.  109 

.078 

.110 

.  103 

.078 

.065 

12 

.068 

.  12b 

.  113 

.  132 

.079 

.079 

13 

.  131 

.  152 

.  153 

.  145 

.132 

.112 

14 

.  128 

.  158 

.  134 

.079 

.068 

.076 

15 

.  127 

.  167 

.  125 

.115 

.049 

.043 

16 

.  162 

.  139 

.  160 

.  142 

.111 

.113 

17 

.  172 

.  151 

.20  8 

.  120 

.077 

.213 

18 

.  123 

.  173 

.  100 

.099 

.146 

.076 

19 

.  164 

.  113 

.096 

.112 

.084 

.107 

20 

.  105 

.  137 

.  108 

.  148 

.  116 

.081 

21 

.  123 

.  178 

.  140 

.  170 

.  130 

.118 

22 

.  192 

.078 

.250 

.  270 

.063 

.05b 

23 

.211 

.  192 

.087 

.  164 

.  113 

.  171 

24 

.244 

.209 

-  195 

.  192 

.  154 

.  148 

25 

.  133 

.  130 

.078 

.  145 

.122 

.133 

26 

.06  5 

.  067 

.  194 

.059 

.056 

.056 

27 

.000 

.000 

.000 

.000 

.000 

.091 

28 

.093 

.  102 

.06  2 

.  115 

.  104 

.083 

29 

.  106 

.  165 

.  137 

.  140 

.113 

.049 

30 

.086 

.  119 

.083 

.  101 

.  120 

.089 

31 

.000 

.000 

.000 

.000 

.000 

.000 

32 

.000 

.000 

.000 

.000 

.065 

.031 

33 

.07  2 

.047 

.044 

.  086 

.043 

.110 

34 

.  146 

.  000 

.20  7 

.074 

.074 

.  148 

35 

.  151 

.  066 

.030 

.  114 

.  029 

.154 

36 

.  15  9 

.  139 

.  139 

.098 

.111 

.094 

37 

.  165 

.  113 

.045 

.  103 

.051 

.076 

38 

.086 

.  102 

.093 

.075 

.073 

.056 

39 

.090 

.  093 

.  123 

.  08  9 

.  1 14 

.094 

Transportation  MOS  ever  all  LCS  for  tne  voluntary  retire¬ 
ments  (type  2  loss)  for  the  grades  of  First  Lieutenant, 
Captain,  Major,  and  Lieutenant  Colonel.  This  illustrates 
the  irregular  distribution  resulting  from  voluntary  attri¬ 
tion.  The  number  of  losses,  even  with  aggregation  ever  LOS, 
is  relatively  small.  The  actual  losses  in  tnis  table  range 
from  zero  to  seven. 


IIBLE  XU 

CENIBAL  ATTBITION  HATES 


AVIATION  (MOS  38) 


GRADE 

197  7 

1978 

1979 

1980 

1981 

1982 

2LT 

-03  0 

.020 

.  031 

.026 

.006 

.017 

ILT 

.05  5 

.082 

.032 

.025 

.029 

.018 

CAPT 

.148 

.  173 

.201 

.157 

.  163 

.096 

HAJ 

.068 

.077 

.  06  9 

.049 

.047 

.045 

LICOL 

.  125 

.  171 

.  144 

.145 

.108 

.  145 

COMBAT 

SUPPORT 

GROUP 

(MOS 

07,13, 

20) 

GRADE 

1977 

1978 

1  979 

1980 

1981 

1982 

2LT 

.00  3 

.025 

.  023 

.009 

.010 

.015 

ILT 

.249 

.240 

.  232 

.187 

.197 

.  142 

CAPT 

.08  2 

.082 

.  104 

.125 

.  106 

.080 

MAJ 

.  123 

.119 

.  107 

.145 

.125 

.080 

ITCOL 

.  136 

.193 

.  114 

.155 

.154 

.114 

GROUND 

COMBAT 

GRCUP 

(MOS  03 

,05, 

10) 

GRADE 

1977 

1978 

197  9 

1980 

1981 

1982 

2LT 

.010 

.016 

.  oil 

.008 

.013 

.011 

ILT 

.23  1 

.190 

.220 

-151 

.  137 

.  121 

CAPT 

.08  1 

.059 

.  074 

.093 

.090 

.065 

MAJ 

.07  4 

.051 

.  060 

.060 

.059 

.032 

LTCOL 

.112 

.  129 

.  131 

.139 

.100 

.  109 

TOTAL  MOS 


GcikDE 

1977 

1978 

1979 

1980 

1981 

1982 

2  IT 

.035 

.032 

.  038 

.025 

.020 

.026 

ILT 

.  17  2 

.161 

.  146 

.118 

.111 

.093 

CAPT 

.  129 

.  128 

.  140 

.  133 

,  123 

.088 

MAJ 

.099 

.095 

.  099 

.034 

.075 

.056 

LICOL 

.  119 

.154 

.  155 

.  146 

.  1 10 

.  130 

A  rurther  disaggregation  tc  include  each  LOS  (not  illus- 
ed)  ,  would  demonstrate  the  effect  of  a  smaller  number  of 
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TABLE  XIT 


HOTOfi  TBANSPOET  (EOS  20)  ?OL.  BET.  (LOSS  TYPE  2) 
CEHTBAL  ATTBITION  BATES 


GRADE 

1977 

1978 

1979 

1980 

1981 

1982 

1LT 

.00  0 

.023 

.00  0 

.0  11 

.000 

.000 

CAPT 

.164 

.049 

.  043 

.047 

.022 

.021 

EAJ 

.000 

.298 

.  169 

.235 

.041 

.  107 

LTCOL 

.000 

.857 

.286 

.222 

.000 

.000 

?.  JAMES-STEm  TECHHIOaES  AMD  TALIDATIOM  PBOCEDOBES 

A.  BACKGBOOHD 

The  first  basic  process  in  statistics  is  the  siuple  act 
of  counting  and  the  second  most  basic  process  is  averaging. 
One  use  of  observed  averages  is  to  predict  or  estimate  unob¬ 
served  guantities.  For  example,  a  baseball  player  who  gets 
eight  hits  in  25  times  at  bat  is  said  to  have  a  .320  batting 
average,  which  is  rhe  ratio  of  8  to  25.  In  computing  rhis 
statistic  an  estimate  is  formed  of  the  player's  true  batting 
anility  in  terms  of  his  observed  average  rate  of  success  at 
the  plate.  A  common  reply  to  the  inguiry  regarding  tnis 
baseball  player's  next  100  times  to  bat  would  be  that  he 
would  probably  get  32  more  hits.  In  traditional  statistical 
theory,  no  other  estimation  rule  is  uniformly  better  than 
the  observed  average.  The  attrition  in  cells  is  modeled  as 
a  set  of  Bernoulli  trials  whose  number  is  the  inventory  of 
that  cell.  The  process  of  choosing  leavers  is  viewed  as 
independent  from  individual  to  individual.  A  convient 
anology  might  be  that  of  a  baseball  players  batting  average 
in  which  each  time  at  bat  is  viewed  as  an  indepenaent 
Bernculli  trial. 

Using  the  above  Bernoulli  (Binomral)  model,  the  rates 
(or  batting  averages)  are  the  Naximum  Likelihood  estimators 
of  the  Bernoulli  parameter.  The  maximum  lixeiihood  proce¬ 
dure  examines  the  likelihood  function  of  the  sample  values 
and  takes  as  the  estimates  of  tne  unknown  parameters  those 
values  that  maximize  the  likelihood  function  [Bef.  4;  p. 
363  ].  These  estimates  are  optimal  (in  particular  admis¬ 
sible)  if  we  are  talking  about  one  cell  (or  one  ball 
player)  .  But  it  is  no  longer  true  if  we  are  dealing  witn 
three  or  more  cells  (cr  ball  players) . 


Charles  Stern  [fief.  5]  showed  in  1955  that  it  is 
possiJble  to  make  a  uniform  improvement  on  the  mazimum  like¬ 
lihood  estimator  (MLE)  in  terms  of  total  squared  error  risk 
when  estimating  several  means  from  independent  normal  popu¬ 
lations.  In  1961*  James  and  Stein  [fief.  6]  presented  to  the 
Berkeley  Symposium  an  estimator  which  shrinks  the  maximum 
likelihood  estimator  toward  the  origin,  and  can  improve 
(i.e.  lower  risk)  on  the  maximum  likelihood  estimator  quite 
substantially  provided  there  are  three  or  more  cell.  To 
continue  the  baseball  example,  if  one  considers  three  or 
more  taseball  players  and  desires  to  predict  future  batting 
averages  for  each  of  them,  then  there  is  a  procedure  that  is 
better  than  simply  extrapolating  from  the  separate  averages. 
Efron  and  Morris  [Ref.  7]  demonstrated  an  application  of  the 
Stein  rule  and  its  generalizations  applied  to  predict  base¬ 
ball  averages,  to  estimate  toxomosis  prevalence  rates,  and 
to  estimate  the  exact  size  of  Pearson's  chi-sguare  test  with 
results  from  a  computer  simulation.  The  result  of  tnese 
applications  was  a  mean  square  error  of  the  Stein  rule  being 
less  than  half  of  that  when  one  uses  the  individual 
empirical  averages. 

The  first  step  in  applying  Stein's  method  is  to  deter¬ 
mine  the  average  of  the  averages.  This  grand  average  or 
grand  mean  must  lie  between  0  and  1.  The  essential  process 
in  Stein's  method  is  the  "shrinking"  of  all  the  indiviaual 
averages  toward  this  grand  mean.  When  the  Stein's  method  is 
applied  to  attrition  rates,  the  results  are  as  follows. 

1.  If  an  aggregated  rate  is  higher  than  the  grand  mean 
then  it  will  be  reduced. 

2.  If  an  aggregated  rate  is  lower  than  the  grand  mean 
then  it  will  be  increased. 

Figure  5. 1  is  a  simple  illustration  of  this  "shrinking" 
concept  utilizing  fictitious  values  for  the  empirical 
values.  In  this  simple  illustration  the  effective  shrinkage 
is  25  percent. 
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JAMES-STEIN  ESTIMATOR  EXAMPLE 


PERCENT  SHRINKAGE 


B 


JAHES-STEIN  ESTIHATOR 


I 


Ihe  James-Stein  estimator  for  each  desired  cell  is  found 
through  the  following  equation. 

1.  Let  'I*  =  the  grand  mean. 

2.  Let  m  =  sample  attrition  rate  for  each  cell. 

3.  Then,  p  =  Y*  ♦  c(m-Y*),  for  a  specialized  value  of 
"c",  is  the  James-Stein  estimator. 

See  Appendix  B,  for  the  James-Stein  estimation  algo¬ 
rithm.  See  Appendix  C,  for  more  mathematical  details  on  the 
equation  used  foe  the  James-Stein  estimator.  The  amount 
(m-Y*)  is  the  amount  by  which  the  sample  attrition  rate 
differs  from  tne  grand  mean  of  the  desired  cells.  The  equa¬ 
tion  thus  states  that  the  James-Stein  estimator  "p*'  differs 
from  tne  grand  mean  by  this  same  quantity  (m-Y’*‘)  multiplied 
by  a  coefficient,  "c".  The  coefficient  "c"  is  the 
"shrinking  factor"  used  in  the  James-Stein  estimation 
process.  If  "c"  was  equal  to  1,  then  the  resultant 

James-Stein  estimator  for  a  given  cell  would  be  identical  to 
the  cell's  sample  attrition  rate.  In  other  words,  p=iu. 
Stein's  theorem  [Hef.  6]  states  that  the  shrinking  factor  is 
always  less  than  1.  Its  actual  value  is  determined  by  tue 
collection  of  all  observed  attrition  rates. 

As  an  example,  if  the  shrinking  factor  "c"  was  equal  to 
.3,  then  each  attrition  rate  will  shrink  about  70  percent  of 
the  distance  to  the  grand  mean.  Three  examples  follow: 

1.  Assume  the  sample  attrition  rate  was  .10,  the  grand 

mean  was  .05,  and  "c"  equal  to  .3.  Then 

p=.  05+ .  3  (.  10- .  05)  .  The  result  is  .065.  By  Stein's 
theorem,  the  attrition  rate  is  best  estimated  by  .065 
rather  than  the  historical  sample  attrition  rate. 
Shrinkage  was  down  toward  the  grand  mean. 

2.  Assume  the  sample  attrition  rate  was  .50,  the  grand 

mean  was  .05,  and  "c"  equal  to  .3.  Then 
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The  esti- 


p=. 05+ .3 (. 50- « 05) .  The  result  is  .185. 
iiiator  is  smaller  than  the  sample  attrition  rate. 
Shrinkage  was  down  toward  the  grand  mean. 

3.  Assume  the  sample  attrition  rate  was  .05,  the  grand 
mean  was  .10,  and  "c"  equal  to  .7.  Then 

p=. 10+. 7 ( . 05- .  10) .  The  result  is  .065.  Shrinkage 
was  up  toward  the  grand  mean. 

There  are  several  expressions  for  the  James-Stein  esti¬ 
mator,  but  all  of  them  have  in  common  a  shrinking  factor 
"c”.  without  exhibiting  any  of  the  various  formulas  for 
"c",  the  following  is  a  description  of  their  general 
behavior.  See  Tanle  XV  for  the  James-Stein  formula  utilized 
in  this  project.  See  Appendix  C  for  further  background  and 
formulas. 

Let  "K"  be  the  number  of  unknown  means;  let  sigma  be  the 
within  population  standard  deviation;  and  let  SSB  be  tne  sum 
of  square  of  cell  means  measured  from  the  grand  mean. 
Holding  sigma  and  “K''  constant  in  the  formula  for  "c",  the 
shrinking  factor  is  affected  by  SSB.  The  shrinking  factor 
"c**  becomes  smaller  as  SSB  gets  smaller.  The  predicted 
means  are  more  affected  by  this  condition  and  the  resultant 
mean  is  closer  to  the  grand  mean.  On  the  other  hand,  as  SSB 
increases,  "c"  increases  and  the  shrinkage  is  less  drastic. 
In  effect  the  James-Stein  procedure  makes  a  preliminary 
estimate  that  ail  the  unobservable  means  are  near  the  grand 
mean.  If  the  data  supports  the  condition  wherein  the  sample 
means  are  not  too  far  from  tne  grand  mean,  then  the  esti¬ 
mates  are  ail  snrunk  further  toward  tne  grand  mean.  If  this 
condition  is  contradicted,  then  the  shrinking  is  minimal. 
The  distribution  of  the  sample  means  around  the  grand  mean 
effects  tne  shrinking  factor.  The  number  of  means  being 
estimated  also  influences  the  shrinking  factor.  If  "K"  is 
rather  large,  the  shrinking  to  the  grand  mean  maybe  more 
drastic.  [Ref.  8:  p.  123] 


54 


■j  « 


[v't’C 


TABLE  XV 
FOBHOLAS 

I.('t  (f  )  inventory  at  the  l.cginning  of  fis.  al  ><  ar  t  (i  I  1' 
Let.  y  riiiniber  of  attritii>n.s  at  any  iiiiie  during  w  ar  i 
l.ci.  It  -  rria.ximum  y  .U.3(.V(/)  •  \(/  -  1)) 

MAXIMUM  LIKELIHOOD 


MIMMAX 


( I  -  \  n 


\  n 


J  AMES-STEIN 

Tlie  following  steps  are  utili/.ed  lo  generate  Jatiies-St ein  altrilion  rates 


S'LEP  1.  Use  a  variance  staLalizing  transform  (Freeman  -  Tukeyj 

[.Sin  ' 


n  •  1 )  (  n  -  1 ) 


D] 


!s'l  LI’  2:  Form  tlie  cell  means,  the  grand  mean  SSH.  and  .sSE 

l.i'l .  i-ell  means 

l.'-i  .\  grand  mean 

I. cl  A  II  n  m  lier  of  I  real  mciil  cells 

S  ri.I’  o  Compute  the  '■ei  of  .lame^-Stein  "siinialor'  in  the  iransforme.i  seal 

,  TiK  :!|  s.'-'i: 

h {T  1 )  2  > ' H 

r  -V  c  ( \  1 

s  i  l  l'  I  Insert  I  Ilf  transform  n.  produce  the  alirition  rate,  ’r" 

r  r  t  .''III  ( - -  I 
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Which  set  of  values,  "m”  or  *'p”,  is  the  better  indicator 
of  the  attrition  rate?  In  order  to  answer  that  question  one 
would  have  to  know  the  true  attrition  rate.  This  true 
average  shall  be  designated  with  "TH”  in  the  main  text  of 
this  paper  ana  by  tne  Greek  letter  theta  in  the  Appendices. 
It  IS  the  probability  that  an  officer  would  leave  any  time 
and  is  actually  an  unknowable  quantity.  Although  "TH”  is 
unooser vable,  an  approximation  can  be  made  utilizing 
historical  data. 

One  method  of  evaluating  the  two  estimates,  i.e.,  the 
James-Stein  estimates  and  the  observed  attrition  rate  m,  is 
ny  simply  counting  their  closeness  to  "TH".  Consider  the 
closer  the  estimator  is  to  “TH''  as  a  success  and  the  farther 
from  "TH"  as  a  failure.  The  number  of  successes  of  one 
estimator  could  then  be  compared  to  the  numoer  of  successes 
of  the  ether  estimator. 

A  more  quantitative  method  of  comparing  tne  two  esti¬ 
mates  is  through  the  total  squared  error  estimation.  inis 
is  measured  by  first  determining  the  actual  error  of  each 
prediction  ,  given  by  (TH-m)  and  (TH-p)  for  each  cell.  ^a.ch 
of  these  quantities  is  then  squared  and  the  squared  values 
are  added  up.  The  observed  averages  “m"  would  have  a  total 
squared  error  denoted  by  Em,  whereas  the  squared  error  of 
the  James-Stein  estimators  would  be  denoted  by  £js. 
Comparing  Em  and  Egs  by  the  ratio  R=Em/Sjs,  then  the 
James-Stein  method  is  “R"  times  as  accurate  as  the  observed 
attrition  rate  "ra". 

Employing  the  ideas  borrowed  from  statistical  decision 
theory,  the  estimators  can  be  compared  through  a  risk  func¬ 
tion.  The  risk  function  of  a  decision  rule  is  the  expected 
loss  (over  the  sample)  incurred  using  the  decision  rule 
[fief-  9:  p.  8].  The  above  techniques  are  similar  to  the 
squared-error  loss  method  [Hef.  9:  p.  9]. 
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Iberetore,  the  risk  function  is  the  expected  value  of 
the  squared  error  for  every  possible  value  of  "IH”.  Stein's 
theorem  is  concerned  with  estimation  of  several  unknown 
means  [fief.  7].  No  relation  between  the  means  need  be 
assumed  and  indeed  are  assumed  to  be  independent  of  one 
another.  The  historical  development  of  statistical  theory 
from  Gauss  through  decision  theory  argues  that  the  average 
is  an  admissible  estimator  as  long  as  there  is  just  one  or 
two  means  to  be  estimated.  An  estimator  for  a  parameter  is 
admissible  if,  according  to  a  specified  criterion,  taere 
exists  no  other  estimator  that  is  better  than  it  for  all 
possible  values  of  the  parameter.  Stein  in  collaboration 
with  James  showed  that  when  the  number  of  means  exceeds  two 
that  estimating  each  of  the  means  by  their  own  average  is  an 
inadmissible  procedure.  No  matter  what  the  values  of  true 
mean,  there  are  estimation  rules  with  smaller  total  risk 
[Ref.  6]. 

The  risk  function  for  the  James-Stein  estimator  is  less 
than  the  risk  function  for  the  sample  means  irrespective 
what  the  true  values  of  the  means  "TH”  happen  to  £e.  The 
reduction  of  risk  can  be  substantial,  particularly  when  the 
number  of  means  is  relatively  large.  Tne  risk  of  the 
James-Stein  estimator  is  smallest  wnen  all  the  true  means 
"TH”  are  equal.  As  the  true  means  increase  in  variation 
from  one  another  the  risk  of  the  estimator  increases, 
approaching  the  value  of  the  observed  averages  but  never 
quite  equaling  it.  The  James-Stein  estimator  does  sunstan- 
tiaiiy  tetter  than  the  averages  only  if  the  true  means  lie 
near  each  other.  The  James-Stein  estimator  does  at  least 
marginally  better  no  matter  what  the  true  means  art  [fief.  8: 
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The  model  for  the  James-Stein  estimator  used  thus  far 
shrinks  the  observed  averages  to  the  grand  mean.  This  is 
not  the  only  possible  procedure.  Other  models  for  the  esti¬ 
mator  dispense  with  the  grand  mean  entirely.  The  observed 
attrition  rates  do  not  depend  on  a  choice  of  origin.  Before 
Stein  discovered  his  method  it  was  generally  accepted  that 
such  “invariant"  estimators  must  be  preferable  to  tnose 
whose  predictions  change  with  each  choice  of  an  origin.  If 
the  origin  or  zero  is  chosen  as  the  grand  mean  then  the 
terms  containing  the  grand  mean  would  be  removed  from  the 
equation.  In  this  case  the  James-Stein  estimator  would  be 
p=c(m).  The  estimation  process  is  now  complicated  by  the 
fact  that  the  shrinking  factor  “c“  would  be  different  tor 
each  cell.  The  shrinking  factor  is  dependent  on  the  stan¬ 
dard  deviation  of  the  sample  attrition  rates.  A  large  stan¬ 
dard  deviation  implies  a  high  degree  of  randomness  or 
uncertainty  in  the  sample  data.  If  the  sample  attrition  is 
large,  it  can  be  attributed  to  random  fluctuations  rather 
than  to  an  actual  large  value  of  the  true  mean  “TH“.  Inus, 
application  of  a  small  shrinking  factor  would  reduce  tne 
value.  [Hef.  7:  pp.  123-4] 

There  is  one  purpose  for  which  the  measured  sample 
attrition  rate  may  well  be  superior  to  the  James-Stein  esti¬ 
mator:  wnen  a  single  cell'  is  considered  in  isolation.  The 
James-Stein  method  gives  better  estimates  for  a  majority  of 
cells  and  it  reduces  the  total  error  of  estimation  for  tne 
sum  of  ail  the  ceils.  Estimating  the  true  mean  for  an 
isolated  cell  by  Stein’s  method  creates  serious  errors  when 
that  mean  has  an  atypical  value.  The  reduction  of  the  risk 
is  more  predominant  in  a  Homogeneous  setting.  Tne  inclusion 
of  an  atypical  or  ncnhomogeneous  mean  in  the  estimation 
process  would  increase  the  risk. 


The  formula  cor  the  James-Stein  is  similar  to  the 
Bayesian  equation  of  Z=m+C(y-ni).  Here  "y”  is  the  sample 
meau/  "m"  is  the  mean  of  the  prior  distribution.  The 
shrinking  factor  “C"  is  different  in  that  it  depends  on  the 
standard  deviation  of  the  prior  distribution.  [Bef.  9:  pp. 
99-110] 

The  James-Stein  procedure,  however,  has  one  impor¬ 
tant  advantage  over  the  pure  Bayesian  method.  The 
James-Stein  estimator  can  be  employed  without  knowledge  of 
the  prior  distribution.  Toe  James-Stein  estimator  can  be 
regarded  as  an  empirical  Bayes  rule.  The  empirical  Bayes 
approach  uses  historical  data  to  estimate  the  prior 
distribution  [Bef.  9:  p.  117]. 

2.  Modifications  to  the  James-Stein  Process 

The  James-Stein  estimation  process  can  be  modified 
using  the  Sfron-Morris  limited  translation  version  of  the 
estimator  [Bef.  10].  This  modification  ensures  that  tne 
estimator  of  "TH"  is  not  shifted  so  far  from  the  sample  mean 
tnat  the  estimator  is  inconsistent  with  the  sample  mean. 
The  modification  is  a  compromise  between  the  James-Stein 
estimator  and  tne  Maximum  Likelihood  Estimatcr  which  has 
good  individual  properties.  The  compromise  would  follow  the 
James-Stein  rule  as  closely  as  possible  subject  to  a  fixed 
constraint  cc  now  far  the  estimator  is  allowed  to  deviate 
from  the  MLE.  Such  "limited  translation  estimators”  were 
discussed  in  a  Bayesian  framework  in  the  Efrcn-Morris 
article  [Bef.  11].  Tney  showed  that  it  is  possible  to 
considerably  reduce  the  maximum  possible  risk  for  any  compo¬ 
nent  wnile  reducing  the  ensemble  savings  of  the  James-Stein 
estimator  by  as  little  as  5  or  10  percent  [Bef.  10]. 


D.  VALICATIOH  TESTS 


A  validatioa  test  was  conducted  to  evaluate  the  effi¬ 
ciency  of  tne  James-Stein  shrinkage  estrinator.  The  test  was 
conducted  as  follows: 

1.  Select  a  grade  within  an  Occupational  Group  to 
examine.  The  resultant  desired  data  array  will  be 
three  dimensional  (years,  LOS,  MOS) . 

2.  Let  "r"  stand  for  LOS,  then  i=0,...,30. 

3.  Let  "j”  stand  for  MOS,  (tnen  values  of  j's  depend  on 
whicn  MOS  group  is  being  analized) . 

4.  Let  Dij  =  Incidence  Matrix  of  nonstructural  zeroes, 
which  is  the  same  for  all  years.  Dij  =  1,  if  cell  is 
member  of  the  feasible  region;  Dij  =  0,  otherwise. 

5.  Let  K  =  number  of  feasible  cells,  i.e.  ,  sum  of  all 
Dij . 

6.  Let  lij  =  number  of  leavers  in  cell  (i,j). 

7.  Let  t  =  1,...,T:  where  T  =  number  of  years  of  data 

used  to  create  the  estimator. 

8.  Let  Nij  =  Inventory  in  cell  (i»j)  = 

Max((N(t)  <-N(t  +  1))/2,I{t)). 

The  validation  procedures  used  t=1...4  to  compute  the  empir¬ 
ical  estimates  and  used  t=5,6  for  validation  purposes. 
Three  estimation  methods  were  employed:  James-Stern, 

Minimax,  ana  Maximum  Likelihood. 

1  Steps : 

The  following  steps  were  utilized  to  prepare  the 
data  for  the  validation  procedures.  These  steps  included 
the  transformation  of  the  given  data  of  leavers  and  computed 
inventory  according  to  the  variance  stabilizing  arc  sine 
transformation  listed  in  Table  XV.  See  Appendix  D,  Figure 
D.18  for  tne  APL  listings  of  the  following  functions. 


cO 


•  let  ISi j  (t)  =  yij(t)  BINPBEP  Nij(t),  where  "BINPREP*' 

prepares  the  Freefflan-Tukey  version  of  the  arc  sine 
transformation  for  binomial  data. 

•  let  T  =  4,  years  of  data  desired. 

•  let  Jij  =  Dij  JAMES  ISij (t) ,  where  "JAMES"  returns  a 
James-Stein  estimator  for  the  means  of  the  cells  in  the 
last  two  dimensicns  of  "IS"  while  being  screened  by  the 
incidence  matrix  "D"  of  nonstructural  zeroes. 

•  let  Pij  =  Jij  BINCONV  Nij(t),  where  "BlhCONV"  inverts 
the  arc  sine  transformation  used  in  preliminary  step 
cne. 

2 .  Validation  Procedures : 

Ihe  following  procedures  were  utilized  to  validate 
the  effectiveness  of  the  the  James-Stein  estimation  process. 

•  let  I5*ij  =  ISij(r*l),  waich  is  assumed  to  be 

distributed  NormaI(Jij, unknown  variance). 

•  let  NOB  =  IS*ij  -  Jij,  where  Jij  was  derived  in 

preliminary  step  three. 

•  Cse  the  distribution  fitting  capaoility  of  GBAFSTAI  or 
any  other  comparable  graphical  display  package,  to 
compare  NOR  to  the  Normal  distribution.  Desired  output 
is  a  comparison  aistogram  of  the  data  to  a  fitted 
Normal  density,  compare  sample  distribution  to  fitted 

'  Normal  distribu tion,  normal  probability  plot  with 

data,  and  possibly  the  survivability  curve  fit.  These 
plots  can  all  be  plotted  on  one  sheet  of  output. 

•  Compute  test  statistics  comparing  data  fitted  to  a 
Normal  distribution. 

•  Compute  the  figure  of  merit  from  the  test  statistics 
using  the  sum  of  the  square  of  the  sample  mean  and  the 
square  of  the  sample  standard  deviation. 


3.  M3.P3.max  Estimator 

The  foUowlpg  procedures  were  utilized  tc  compare 
the  James-Stein  estimator  to  the  Minimax  estimator.  See 
Table  XV  for  the  formula  for  the  Minimax  estimator.  An 
simple  example  of  the  possible  values  resulting  from  the  use 
of  this  formula  is  provided  in  Table  XVI,  to  illustrate  the 
shrinkage  characteristics  inherent  in  its  use  for  the 
purpose  of  projecting  rates  which  have  extreme  empirical 
values. 


TABLE  IVI 

HIMIHAX  EXAMPLE 

LEAVERS  INVEHTORY 

VALUES 

MINI.MAX 

ESIIMATE 

0 

0 

.  50 

0 

1 

.25 

1 

1 

.75 

0 

4 

.  17 

1 

4 

.  33 

See  Appendix  D,  Figure  D. 19  for  the  APL  listings  of  tnis 
function. 

•  let  PMMij  =  Yij(t)  MINMAX  Nij  (t)  ,  where  «MIN:1AX'' 
returns  the  MIMMAX  estimates  for  the  binomial. 

•  let  ISiS  =  Nij  ABCSIN  PMMij,  where  "ARCSIN”  returns  the 
inversfc  sine  transformation  for  use  when  the  success 
probability  PMMij  is  estimated  directly  (e.g.,  by 
MI  NMAX)  . 

•  let  NOlMM  =  IS*  -  ISM,  where  IS*  is  from  validation 
procedures  step  one. 


•  Us€  the  distribution  fitting  capability  of  GRAFSTAI  or 
any  other  comparable  graphical  display  package,  to 
compare  NOBHN  to  the  Normal  distribution  similar  to  the 
validation  procedure  three. 

•  Compute  test  statistics  comparing  data  fitted  to  a 
Normal  distribution. 

•  Compute  the  figure  of  merit  from  the  test  statistics 
using  the  sum  of  the  square  of  the  sample  mean  and  the 
square  of  the  sample  standard  deviation. 

•  Maximum  Like  lihood 

The  natural  empirical  estimator  without  shrinkage  is 
the  Maximum  Likelihood  estimator.  Comparison  of  the 

James-Stein  process  to  the  Maximum  Likelihood  process  will 
illustrate  the  efficiency  of  tne  James-Stein  estimation 
process. 

•  let  Y  =  sum  of  Yij(t)  for  t=l  to  T.  Where  T  =  4. 

•  let  N  =  sum  of  Nij(t)  for  t=1  to  T.  Where  T  =  4. 

•  let  ISAij  =  Y  BINPREP  N. 

•  let  NORA  =  IS*  -  ISAij,  where  IS*  is  from  the  valida¬ 
tion  procedures  step  one. 

•  Ose  the  distribution  fitting  capanility  of  GRAFSIAT  or 
any  other  comparaoie  graphical  display  package,  to 
compare  NORA  to  the  Normal  distribution  similar  to  the 
validation  procedure  three. 

•  Compute  test  statistics  comparing  data  fitted  to  a 
Normal  distribution. 

•  Compute  the  figure  of  merit  from  the  test  statistics 
using  the  sum  of  the  square  of  the  sample  mean  and  the 
square  of  the  sample  standard  deviation. 


E.  OITA  AHALYSIS 

A  data  analysis  was  conducted  utilizing  the  three 
comparison  mentioned  above. 

1 .  The  James-St ein  Estimator 

Four  years  of  data  (1977-1930)  was  used  to  construct 
the  James-Stein  estimator.  The  leavers  and  inventory  were 
transformed  with  the  Freeman-TuJcey  transformation.  The 
projected  year*s  (1981-1903)  leavers  and  inventory  were 
constructed  and  also  transformed  with  the  Freeman-Tukey 
transformation.  The  difference  between  the  James-Stein 
estimator  and  the  projected  years  data  was  plotted  to 
compare  the  fit  to  a  Normal  distribution.  The  figure  of 
merit  (FOd)  was  computed  from  the  inputs  taken  from  the 
statistical  table. 

FOM=  (squared  mean)  (variance)  . 

Table  XVII  illustrates  the  results  of  the  computa¬ 
tion  cf  the  James-Stein,  Hinimax,  and  Maximum  Likelxhood 
FOM*s.  The  data  organization  format  of  the  table  is 
provided  for  ease  of  interpretation  of  the  table.  The  table 
demonstates  the  calculation  of  the  FOM  for  the  years  1981, 
1982,  and  1983.  The  separate  cases  evaluated  are  identified 
by  rows  which  represent  the  grades  of  First  Lieutenant  and 
Lieutenant  Colonel  in  each  of  the  four  dOS  groups  utilized 
in  this  project. 

2  .  Ih e  Minima  Estimator 

The  leavers  and  inventory  of  the  four  years 
(1977-1980)  were  used  to  construct  a  minimax  estimator.  See 
Table  XV  for  the  formula  utilized  to  derive  this  estimator. 
The  result  was  compared  to  the  fit  of  the  Normal  distribu¬ 


tion  in  a  similar  manner  as  discussed  above. 


The  FOM  was 


computed  as  shown  above  and  further  illustrated  in  Table 
XVII. 

3 .  Maximum  Likelihood 

The  last  comparison  was  using  the  leavers  and  inven¬ 
tory  of  the  four  years  (1977-1980)  and  only  using  the 
Freeman-Tukey  transf ocmation.  The  result  is  the  Maximum 
Likelehood  estimator  in  the  transformed  space.  The  algo¬ 
rithm  is  the  same  as  James-Stein  estimator  described  above, 
only  the  James-Stein  calculation  is  by-passed.  The 
projected  years  data  was  used  to  compare  the  resultant 
difference  with  the  fit  to  the  Normal.  The  FOH  was  computed 
as  shewn  above  and  further  illustrated  in  Table  XVII. 

4.  Figure  of  Merit 

Ihe  figure  of  merit  is  representative  of  the  effi¬ 
ciency  of  the  estimator.  Since  the  FOM  is  computed  in  the 
transformed  and  Normalized  scale,  the  smaller  the  FOM,  the 
oetter  the  estimate.  The  ratios  illustrated  in  columns  3 
and  5  of  Table  XVII  illustrate  (in  most  cases)  how  many 
times  larger  the  other  estimating  technigue's  FOM  is  when 
compared  to  the  James-Stein  estimator  technique.  One  can 
note  that  in  one  case  the  use  of  the  James-Stein  estimator 
technique  resulted  in  an  improvement  of  as  much  as  30  times 
over  the  Maximum  Likelihood  process  and  22  times  over  the 
Minimax  process.  These  results  are  typical  of  other  cases 
executed  by  the  author,  but  not  shown  for  redundancy  sake. 

F.  ALDITIONAL  COBPABISOBS 

The  above  risk  ccoparisons  were  all  performed  in  the 
transformed  (Freeman-Tukey)  space.  Since  this  distorts  the 
scale  it  seems  wise  to  perform  some  additional  comparisons 
(sum  of  squares  error)  in  the  original  (untransformed) 


TABLE  XVII 

F160BE  OF  HEBIT  COHPAfilSONS 


♦♦♦♦♦♦♦*♦♦♦  DATA  ORGANIZATION  FORMAT  ♦♦♦♦♦♦♦♦♦♦♦ 
♦♦♦  COLOMN  ♦♦♦  *♦*  REMARK  ♦♦♦ 


JAMES 
MAXLIKE 
HAXIIKE/JAHES 
MINIMAX 
MINIMA X/JAMES 

♦  ♦♦  ROW  ♦♦  = 


**** 


FOM  Of  data  using  James- Stein 
FON  of  Maximum  Likelihood 
Ratio  of  Haxlike  divided  by  James 
FOM  of  data  using  Minimax 
Ratio  of  Minimax  divided  by  James 

♦♦♦  REMARK  ♦♦♦ 

1ST  LT  in  Aviation 

1ST  IT  in  Combat  Service  Support 

1ST  LT  in  Combat  Support 

1ST  LT  in  Ground  Combat 

LTCOL  in  Aviation 

LTCOL  in  Combat  Service  Support 

LTCOI  in  Combat  Support 

LTCOL  in  Ground  Combat 


**** 

THESE 

ARE  THE  FIGURES  OF  HER 
MAX LIKE/ 

IT  FOR  1981  ♦♦♦ 

MINIMA 

RCW 

JAMES 

MAXLIKE 

JAMES 

MINIMAX 

JAMES 

1 

1.69 

51.47 

30.46 

37.68 

22.30 

z 

0.71 

1.60 

2.25 

0.78 

1.10 

3 

1.20 

3.69 

3.08 

2. 10 

1.75 

4 

1.  19 

5.47 

4.59 

3.81 

3.20 

c 

w 

3.05 

13.64 

4.47 

8.92 

2.92 

6 

0.22 

0.42 

1.91 

0.23 

1.05 

7 

0.31 

2.34 

7.55 

1.16 

3.74 

8 

0.79 

6.81 

8.62 

3.73 

4.72 

THESE 

JAMES 

4.37 

0.75 

1.53 

1.89 

5.80 

0.23 


ARE  THE  FIGURES  OF  MERIT  FOR  1  982 

MAXLIKE/  MINIMAX/ 

MAXLIKE  JAMES  MINIMAX  JAMES 


47.05 
1.55 
3.71 
5.84 
17.  28 
0.54 


10.77 
2.07 
2.  42 
3.09 
2.  98 
2.35 


34.29 

0.76 

2.19 

4.21 

13.08 

0.33 


7.85 

1.01 

1.43 

2.23 

2.2b 

1.43 


7 

0.59 

2.  60 

4.41 

1.43 

2.42 

8 

1.13 

8.  34 

7.38 

4.92 

4.35 

**** 

THESE 

ARE  THE  FIGURES  OF  MERIT  FOR  1  983  **♦ 

MAXLIKE/  MINIMA 

RON 

JAMES 

HAXLIKE 

JAMBS 

MINIMAX 

JAMES 

1 

4.64 

43.  16 

9.28 

30.65 

6.63 

2 

0.91 

1.64 

1.80 

0.86 

0.95 

3 

1.  17 

3.67 

3.  14 

2.02 

1.73 

4 

2.85 

5.67 

1.99 

4.  21 

1.48 

c 

5.49 

11.93 

2.  17 

8.44 

1.54 

6 

0.2b 

0.56 

2.  15 

0.34 

1.31 

7 

0.62 

2.76 

4.45 

1.53 

2.47 

8 

1.57 

8.55 

5.  45 

5.16 

3.29 

scales.  Accordingly  the  deviations  of  the  empirical  attri¬ 
tion  rates  in  the  validation  time  frame  from  the  estimates 
(both  James-Stein  and  Narimam  Likelihood)  are  computed  and 
compared.  Specifics  follow: 

The  preliminary  steps  are  as  follows: 

•  Actual  rates:  empirical  attrition  rates  are  the  central 
attrition  rates  for  the  projected  period. 

•  James-Stein  Projected  rates:  projected  attrition  rates 
utilizing  the  James-Stein  estimation  process  previously 
defined. 

•  Let,  (1L£=  Maximum  Likelihood  Estimates,  the  historical 
empirical  attrition  rates  are  assumed  to  remain  the 
same  throughout  the  projected  period. 

•  let,  SEP=  The  squared  differences  between  the 
James-Stein  projected  attrition  rates  and  the  actual 
rates. 

•  let,  SEM=  The  squared  differences  between  the  MLE  and 
the  actual  rates. 

•  let,  SSE(MLE)=  The  sum  of  squared  differences  of  SEM. 

•  let,  SSE  (P)  =  The  sum  of  squared  differences  of  SEP. 

•  let,  SEP<SEM  be  the  number  of  cells  for  which  SEP  was 
less  than  SEN.  (For  convenience  of  column  heading, 
inequalities  are  allowed  to  represent  numbers.) 

•  Let,  SEM<SEP  be  the  number  of  cells  for  which  SEJ  was 
less  than  SEP. 

•  P.OS  Groups  evaluated  were:  CS  -  Combat  Support  Group; 
GC  -  Ground  Combat 

•  Banks  evaluated  were:  ii  -  First  Lieutenant;  ITCOL  - 
lieutenant  Colonel 

1 .  Subtest  1. 

First  it  was  assumed  that  all  cells  in  the  feasinle 
region  were  eligible  fcr  forecasting  the  projected  attrition 
rates.  The  results  of  subtest  1  are  illustrated  in  Table 


XVIII.  The  sum  of  squared  differences  for  the  HLE  and  the 
James-Stein  projected  ’rates  for  the  first  and  last  groups 
are  competitive.  One  has  to  keep  in  mind  that  the 
Jame-Stein  estimation  process  will  result  in  a  value  greater 


than  zero  given  a  cell  with  an  original  value  of  zero. 
Therefore,  when  all  the  cells  in.  the  feasible  region  are 
considered  eligible  fcr  projection,  the  empirical  values  of 
zero  result  in  a  value  closer  to  the  grand  mean.  The  value 
of  the  sums  of  the  differences  between  the  James-Stein  esti¬ 
mation  process  and  the  actual  rates  can  be  understandably 
larger  than  the  sums  of  the  differences  between  the  MLS  and 
the  actual  rates  when  there  are  numerous  zero  cells 
originally  in  the  feasible  region. 

I - - - 1 

TABLE  XVIXl 
SOBTEST  1 


Feasible 


Group  Bank 

Cells 

SEP<SEH 

SEMCSEP 

USE  (aiE) 

SSE(?) 

C£ 

LT 

57 

15 

42 

1.52 

1.91 

GC 

LI 

57 

22 

35 

0.33 

3.99 

C£ 

LTCOL 

48 

12 

36 

1.29 

2.93 

GC 

LTCOL 

46 

21 

27 

2.05 

0.66 

2 .  Subtest  2 

A  second  comparison  assumed  a  conditioning  that 
considered  only  the  cells  in  which  there  was  an  actual 
attrition  rate  greater  than  zero.  In  other  words,  what  was 
the  projected  attrition  rate  given  an  actual  less  occured? 
This  created  a  subfeasiole  region  whicn  was  smaller  than  the 


whole  feasible  region.  The  results  of  subtest  2  are  illus¬ 
trated  in  Taole  XIX.  The  sum  cf  squared  differences  for  the 
James-Stein  projected  rates  in  this  subtest  are  much  lower 
than  the  HLE. 


T&BIE  XIX 
SOBTEST  2 


Sur- 

Feas icle 


rcup  Rank 

Cells 

SEPkSSM 

SEM<SEP 

SSE  (MLE) 

SSE(P)  i 

j 

C£ 

LT 

26 

15 

1  1 

1.52 

0.27 

GC 

LT 

29 

22 

7 

0.33 

0.  10  1 

CS 

LTCOL 

16 

12 

4 

1.29 

0.  14  i 

GC 

LTCOL 

28 

21 

7 

2.05 

0.12  1 

i 


VI.  CONC 10 SIGNS  AND  BECOMMENDATIONS 


A.  SOHHAfiX 


The  purpose  of  this  study  was  to  demonstrate  the  appli¬ 
cation  of  the  James-Stein  and  other  shrinkage  type  parameter 
estimation  schemes  for  the  goal  of  generating  manpower  loss 
rates  witnin  the  USMC  officer  force  structure.  This  thesis 
contains  comparisons  of  performance  of  James-Stein,  Minimax, 
and  Maximum  Likelihood  estimations  of  Marine  Corps  officer 
attrition  rates  for  select  cell  aggregates.  The  very  large 
number  of  cells  within  the  USMC  officer  force  structure 
leads  tc  the  condition  that  empirical  attrition  rates  are 
unstable.  This  problem  is  compounded  by  the  fact  that  the 
cell  probabilities  are  small.  Most  rates  are  less  than  10^ 
and  virtualxy  none  are  greater  than  20%  at  disaggregated 
levels.  Further  difficulties  are  present  because  some  of 
the  inventory  cells  are  empty  for  structural  reason  while 
otuers  are  empty  by  chance.  Therefore,  the  small  inventory 
cells  draw  especial  attention.  It  has  been  illustrated 
witnin  tnis  project  tnat  improvement  can  be  attained  by 
application  of  the  James-Stein  and  Minimax  shrinkage  methods 
ratner  than  the  more  natural  Maximum  Likelihood  estimation 
process.  It  is  important  to  note  tnat  the  James-Stein  ana 
Minimax  snrinxage  schemes  seem  to  compete  for  appropriate 
attrition  rate  generation. 


B.  CCNCIUSIONS 


The 

usef  u  1 
ment  sh 
in  man 
exploit 


shrinkage  schemes  employed  herein  offer  powerful  and 
methods  for  generating  attrition  rates.  This  employ- 
juld  lead  to  much  lower  costs  that  accrue  from  errors 
power  planning.  However,  tne  particular  way  to 
them  is  yet  tc  oe  determined. 
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The  application  of  tie  James-Stein  shrinJeage  technique 
has  exposed  a  number  of  problem  areas  which  could  be 
expanded. 

Ibe  aggregation  problem.  Aggregation  of  cells  (cross 
classified  by  grade.  Military  Occupational  Specialty,  length 
of  service,  etc)  into  sets  whose  attrition  behavior  is  homo¬ 
geneous  and  of  low  internal  variability  should  lead  to  the 
most  useful  overall  performance.  One  method,  not  attempted 
by  the  author,  is  to  aggregate  over  MOS  types  in  which  the 
grade  structure  is  siiilar.  For  example,  some  MCS's  are 
restricted  to  warrant  officers  only,  warrant  officers  ana 
limited  duty  officers  only,  unrestricted  regular/reserve 
officers  only,  etc.  Another  method  would  be  to  aggregate 
over  certain  length  of  service  periods.  For  example,  1-4 
years  of  service,  5-10  years  of  service,  10-20  years  of 
service,  and  21-30  years  of  service.  Thus  examining  the 
attrition  rates  in  relationship  to  breakpoints  and  longevity 
in  the  officer's  career  pattern. 

The  limited  translation  problem.  The  shrinkage  estima¬ 
tors  translate  tne  rates  toward  the  aggregate  mean.  Such 
translations  can  appear  to  be  excessive  in  instances  of 
extreme  rates.  This  suggests  that  techniques  that  limit  the 
translation  could  be  profitably  applied.  Current  literature 
contains  limited  translation  methods  wnich  could  be 
examined. 

The  yearly  update  problem.  The  estimated  attrition 
rates  are  based  upon  data  from  several  recent  years.  As 
each  year  produces  a  new  set  of  experiences  there  is  need 
for  a  policy  to  include  new  data  in  the  estimation  scheme 
and  phase  out  the  influence  of  the  distant  past. 

The  validation  prooiem.  The  smallness  of  the  rates 
(i.e.,  0-205?)  places  the  current  project  in  a  range  which 
nas  net  yet  been  treated  successfully  in  tne  literature. 
For  this  reason  it  is  especially  important  to  apply  cross 
validation  procedures  to  score  the  erficacy  of  the  methods. 


Ihe  oultiiiofflial  estimaticn  problem.  The  number  ox 
attritions  in  a  cell  are  placed  into  a  number  o£  (currently 
eight)  disjoint  categories.  This  produces  a  number  o£  small 
sample  size  multinomial  probability  estimation  problems. 
Methodologies  £or  managing  such  problems  have  attracted 
attention  in  the  recent  literature.  Their  use£ulness  could 
be  £urther  examined  in  regard  to  the  officer  attrition  rate 
estimaticn  process 

C.  BiCOHHENDATIOHS 

It  is  recommended  that  further  studies  pursue  the  prob¬ 
lems  identified  above.  It  is  suggested  that  a  new  data  tape 
be  provided  to  the  Naval  Postgraduate  School  which  allows 
the  distinction  between  restricted  and  unrestricted  offi¬ 
cers.  Additionally,  the  data  for  fiscal  years  1984  and  1985 
needs  to  be  provided  in  the  near  future.  In  order  to 
provide  continued  service,  a  regularly  scheduled  yearly 
updated  data  set  should  be  forwarded  to  tne  Naval 
Postgraduate  School  when  available. 


APPESDIX  A 

CENIBAL  ATTBITION  BATE  GBAPHICS 

This  appendix  contains  graphical  illustrations  or  the 
central  attrition  rate  for  USMC  officers. 

Figures  A. 1  to  A. 10  contain  the  central  attrition  rates 
for  each  Military  Occupational  Specialty  distributed  over 
the  length  of  service  period  from  0  to  30.  These  figures 
illustrate  the  aggregate  central  attrition  rate  for  all 
eight  types  of  losses  over  the  six  years  from  1977  to  1982. 
These  figures  demonstrate  the  breakpoints  within  each  MOS 
over  the  distribution  of  length  of  service. 

Figures  A. 11  to  A. 20  contain  the  central  attrition  rates 
for  each  Military  Occupational  Specialty  distributed  over 
the  grades  from  Warrant  Officer  to  Colonel  (labelled  0...9 
on  the  abscissa).  These  figures  also  illustrate  the  aggre¬ 
gate  central  attrition  rate  for  all  eight  types  of  losses 
over  the  same  six  years. 

Figures  A. 21  to  A. 28  contain  the  central  attrition  rates 
for  each  Military  Occupational  Specialty  distributed  over 
the  years  from  1977  to  1982.  Each  figure  nas  five  MOS 
graphs  and  are  dis tinguisharle  by  the  type  of  line  and 
symbol  points  for  each  separate  MOS.  These  figures  illus¬ 
trate  the  aggregate  central  attrition  rate  for  all  eight 
types  of  losses  over  the  six  years  from  1977  to  1982. 

Figure  A. 29  contains  the  central  attrition  rates  for  the 
Aviation  Group  for  the  grades  from  Second  Lieutenant  through 
Lieutenant  Colonel  for  the  years  1977  to  1982.  Figure  A. 30 
contains  the  central  attrition  rates  for  the  Comnat  Support 
Group  for  the  grades  from  Second  Lieutenant  to  Lieutenant 
Colonel  for  tne  years  1977  to  1962.  Figure  A. 31  contains 
the  central  attrition  rates  for  the  Ground  Combat  Group  for 


the  grades  from  Second  Lieutenant  to  Lieutenant  Colonel  for 
the  years  1977  to  1982.  Figure  6T0TM0S  contains  the  central 
attrition  rates  for  the  Total  MOS's  for  the  grades  from 
Second  Lieutenant  to  Lieutenant  Colonel  for  the  years  1977 
to  1982. 


MCIH  or  SERVICE  LENGTH  Of  SERVICE 


Attrition  bates  ny  LOS  for  MOS 


Attrition  Bates  by  LOS  for  MOS*s  09  to  12 


Figure  A. 4  Attrition  Bates  by  LOS  for  HOS* 
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trition  Sates  by  LOS  for  I105*s  17  to  20 


Attrition  Bates  by  LOS  for  MOS 


trition  Rates  by  LOS  for  MOS's  25  to  28 


Attrition  Rates  by  lOS  for  MOS's  29  to  33 


Figure  A. 9  Attrition  Bates  by  LOS  for  MOS*s  34  to  37 
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Figure  £.10  Attrition  Rates  by  LOS  for  UOS 
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.16  Attrition  Bates  by  Grade  for  MOS* 


ttzitiou  Rates  by  Grade  for  HOS’s  25  to  28 
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igure  A. 20  Attrition  Bates  by  Grade  for  noS's  38  to  39 
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tritioD  Rates  by  Year  for  llOS's  26  to  30 


Fig are  A. 27  Attrition  Sates 
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Figure  A. 28  Attritioo  Hates  bj  Year  for  HOS's  36  to  39 
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Figure  i. 30  Ceotral  Attrition  Rates 
for  the  Coabat  Support  Group. 


Figace  A. 31  Centcal  Attrition  Rates 
for  the  Groand  Coabat  Group. 
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Pigace  i.32  Central  Attrition  sates  for  the  Total  dOS's 


JABES- STEIN  ESTIBATOS  ALGOBIIHH 


This  appendix  contains  the  detailed  algorithm  for  the 
James-Stein  estimator  of  attrition  rates  for  the  (ISMC 
officer  force  structure.  The  information  provided  includes: 

•  Assumptions. 

•  Notation. 

•  Steps  to  generate  the  James-Stein  attrition  rates. 


Algorithm;  James-Stein  Estimator  of  Attrition  Rates.  USMC  Officers. 


Segmentation:  Hold  each  grade  fixed  and  aggregate  the  MOS  variable  into  four  main  cells. 
The  four  main  cells  are  as  follows: 

1.  Aviation 

2.  Ground  Combat 

3.  Combat  Support 

4.  Combat  Service  Support  (all  others) 

Notation: 

Let  J  =  the  number  of  MOS  cells  in  the  chosen  aggregate. 

Let  I  =  the  number  of  LOS  cells  in  the  chosen  aggregate  (usually  31). 

Because  the  data  is  a  combination  of  two  types  of  data  (central  and  transition)  and  because 
the  inventory  contains  both  structural  and  sampling  zeroes,  the  following  policies  were  adopted: 

1.  Let,  INV,j{t)  =  inventory  with  LOS  =  i  and  MOS  =  j  at  the  beginning  of  vear  t 

(t  =  l . T) 

2.  Let.  y,y  (t )  =  number  of  attritions  in  cell  (ij)  at  any  time  during  year  t. 

3.  Let,  n,j  (t )  =  maximum  y,y  (t  )  ,  0.5  (/-W'.y  (t )  +  (t  +1))  j 

The  incidence  matrix  D  identifies  the  cells  with  non  zero  inventory. 

Let,  D,j  -  0  if  n,j  (< )  =  0  for  all  t  =  l,...,T 
Let,  D,j  =  1  otherwise. 


The  following  steps  are  utilized  to  generate  James-Stein  attrition  rates. 


STEP  1:  Use  a  variance  stabalizing  transform  (Freeman  -  Tukey). 


=  n.i(t)  - 


(0.5)  (  Sin' 


(2y.,(0) 

(”.,  (0  -  1) 


-  (  Sin  “ '  ; 


2(y.,(0  -  1) 

(%(0  -  1) 


STEP  2:  Form  the  cell  means  and  the  grand  mean. 

For  all  (i.j)  such  that  D„  =  1 

■V  - 

K  -  EE  A; 

I  =  1  ;  -  1 


STEP  3:  Form  SSE,  sum  of  squares  error  and  SSB.  sum  of  squares  between  b\  subtracting 
from  SST,  the  total  sum  of  squares. 

557  SSE  -V.,(t)  X  -D„ 

I  =  I; = II  =  1 
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ssE  =  EEEi^olO  - 
SSB  =  SST  -  SSE 


STEP  4:  Compute  the  set  of  James-Stein  estimators  in  the  transformed  scale. 
T{K-i)  SSE 
^  K{T-l)+2  SSB 


if  £),  ,  =  1 


(X  +  C(Xii-X) 

^  ~  \tnde fined  otherwiee 

STEP  5:  Invert  the  transform  to  produce  the  attrition  rates  r,y . 

[-ifl  *  sin{—  *  0.5) 


^  P 

if  -— <  — 

2  n,  j 


i  'JL 
2^2 

2  '  n.,  "  2 
2  '  n,,  "  2 


■f  ^  ^ 

if  -  — ^ 


..  TT  ^  P 
if  — $ 


«.i  = 

^  (=1 
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LPRMSll  C 
FORMULAS  AMD  PROOFS 


lAis  appendix  contains  various  formulas  and  proofs  for 
the  Janes-Stein  process  utilized  throughout  this  research 
project.  It  is  divided  into  eight  sections  as  follows: 

A.  Simple  Stochastic  Model  for  Central  Data. 

B.  Supporting  Distributions. 

C.  Distribution 

D.  Non  Central  Chi  Square 

£.  Distribution  of  "H"  and  Expectation  of  1/H 

F.  Loss  and  Risk  Functions 

G.  Evaluation  of  E(0/R) 

H.  Inverse  Sine  Transform 

I.  Monthly  Verses  yearly  Rates 


no 


A.  SIMPLE  STOCHASTIC  MODEL  FOR  CENTRAL  DA  TA 


This  section  contains  the  nomenclature  and  characteristics  of  a  simple  stocheistic  model  for 
central  data.  It  deals  with  the  joint  distribution  of  cell  inventor}'  and  number  of  attritions. 


The  following  nomenclature  will  be  defined  in  this  section  as  follows: 
N  =  Inventory  during  period,  which  is  distributed  Poitson  (  A  )  . 
Y  =  Number  of  losses. 


m  =  Central  attrition  rate  =  —  . 

N 


The  conditional  probability  of  losses  given  inventory  is: 


/(J'l  ^  =  ")=  (  ;  )  eMi  -  ©  ) 


y  =  0,...,n 


The  Joint  probability  function  is: 


/  (j/.n)  =  /  (y  '  «)?(«)=  (  J  )  eM  1  -  ©  ) 


0  $  y  ^  n  $  oc 


h{y).P  Y  =  y  =  g  /(„,„)=  g  - 1 - ^(l_e) 

.  =  .  y:  .  =  .  (  n  -  y  )! 


«  -  »  i  ■ 


A  e  ^  *  :(  1  -  e  )  A  ” 


(  n  -  y  )! 


Let  j  =  n  -  y  . 
.  (  A  0)*  . 


1  -  0  )  A  * 


_  (  A  e)»  ,  e 

■  y! 


Therefore  ,  Y  is  Poiston  (  A  0  )  .  which  is  the  marginal  distribution. 

Now.  the  estimated  value  for  the  central  attrition  rate  is: 


.  E  m  --  E  ^  ^  .V)  -  .V)  -  f  A  (-)  -  0 


The  conditional  probability  of  inventory  given  losses  is: 
P(A  =n  V  -  y)  .  <,(n  y]  -  ^ 

n  (y  ) 


(:) 


©Ml-  0 


«  V  (  ©  A  )• 


0  ■-  y  •.  ti  3C 


rW»*»  .?s 


e 


M  i  -  e  ) 


n 


=  y,y-lv. 


X 


A 


This  is  a  shifted  Poisson,  therefore 

N  -  y  is  Poitton  (  A  (  1  0  )  )  . 

Characteristics  of  N  -  y  are: 

MEAN:  E\N  \  Y  =  y]  =  E\N  -  y  +  y  \  Y  =  y] 
=  E\N  -  y\  Y  =  y\-i-  E\ 

-  A  (  1  -  0  )  +  V 

VARIANCE  :  Var  \N  \  Y  =  y\=  Var  [AT  -  y  j  T  = 

COVARIANCE:  Cov  {Y,N)  =  E'Y  N'  -  E  Y 

-  E\E{Y  N  \  N]  - 
=  £!iV  E{Y  ;  N)] 
=  E[N{N  ©  )]  - 
=  £W*©  -  © 

=  0  (  £  (.V*)  -  A 
=  0  (Vor(Ar)) 

=  0  A 


CORRELA  TION ;  Cor 


Cov(Y.N] 

SD{Y)SD{N) 

0  A 

v  0  A  \  ’  A 


-  V"  0 


Therefore. 

\  -  Inventory  in  a  cell  is  distributed  Pouton  (  A 

Y  -  Losses  in  a  cell  is  distributed  Poisson  (  ©  A 

V  N  is  distributed  Binomial  (  ,V .  0  )  . 

0  =  Central  attrition  rale. 

N  Y  is  distributed  Y  ~  Poisson  (  A  (  J  -  ©  ) 


y\y  =  y] 

V]=  A  (  1  -  ©  ) 

; 

(0  A  )  (  A  ) 

-  ©  A  * 

9  A  * 

A^ 


I 

) 


B.  SUPPORTING  DISTRIBUTIONS. 

The  following  sections  of  this  nppendix  cont&in  a  series  of  proofs  to  the  following  statements 
based  on  the  assumption:  Y  is  distributed  Normal  (0,7)  .  where  I  is  the  Identity  matrix. 


i)  U  = 


is  distributed  Normal  ( :  !  6;  ;  ,1)  . 


A'  V  A  A  AA^ 

2  ^ -  I  |e|  |»  =  ^"^1  lei  [■  “ (0-^ -  I  |e|  ' 

See  section  C  for  the  following  statements: 
iii)  V  =  I  \  Z  \  I  ’  is  distributed  At(x-i)  • 
iv)  W  =  £/*+  V  =  I  i  Y\  I  *  . 

PROOFS. 

i)  U  is  independent  of  all  components  of  Z.  In  other  words  all  Cov  {U  ,Z)  -  0  .  Proof  is 
as  follows. 

Cov{U,Z)=  Cov(U,Y-  I  I  I  )  =  Cov{U,Y)-  ,  -  Cov(U,U) 


Now.  Coe(t^  Y)  =  = 


E©,  Coe(Y,  Y)  = 


=  0  .  Q.E.D. 


Therefore,  Cov(U,Z)  =  f ,  , - ^ -  =  0  .  Q.E.D. 

I  i  I  I  0|  i 

ii)  The  distribution  for  Z  is  Normal  and  the  parameters  are  shown  as  follows. 


1.  ft,  2 

,.l  I  ■  0l  i 


E(Z)  =  £(Y)-  -  ^~£(Y)  =  ©- 


©©'© 


E(ZZ  ')-(EZ){EZ)'  =  E{ZZ  ')-0  =  E{ZZ  ')  =  Cov{Z,Z) 


Cov  (Z  ,Z  )  -  Cov  (  Y  -  ■ 


I  I  ©'  ' 

'  !  I  I 


=  Cov  (  Y  ,  Y )  -  Cov  (- 


.Y)-Coe(  Y.- 


Now  Cov  (  1' .  }■ )  =  /  ,  where  I  is  the  identity  matrix  since  ^  is  distributed  .Normal  (©.7  ) 
Also. 

..  ,  re  .  re  re  ,,  re  ,  re 

Co.l  ^  .  -I  El—  H  ^  I  H-  )  £(— g- 


=  £(£='- 


=  EUP)- 


=  (1^  © 


©0' 

;  ©I  2 


KM 


=  ©©'{— 


6 


e 


i) 


=  ee'(- 


’  ©, 


1  -i) 


0©' 


UQ  UQ 

Since  Cov  (  K ,  ,  ,  ^ ,  ■ )  and  Cov  (  ^  ^  ,  y )  are  transposes  of  each  other  and  covari- 


© 


I  I  i  I  i  I  I 

ance  matrices  are  symmetric,  the  result  is  as  follows  using  a  previous  result  for  Cov  [U  ,Y). 


.,Y)-  -2  g— Coe(ir.y)=  -2— 


©I 


I  I  ei 


©I 


Therefore, 


Cov{Z,Z)  =  1-2- 


©©' 


©©' 


and  Z  is  distributed;  Normal  (0,7 


©I  ' *  i  !  ©I  1  * 

©© 


=  7 


©©' 


:  ©i  !  * 


i  ©1 


■) 


%] 


C.  DISTRIBUTION  "  V” 

This  section  is  an  continuation  of  section  B  and  shows  V  -  Z  '  Z  Z  '  ^is  distri- 

K-l 

buted  X(jr-i)  =  2  Normal  (0,1)*- 
1 

PROOFS. 

iii)  First,  construct  an  orthonormal  transformation  P,  where  P  is  K  by  K  dimensional  and 
P  'P  =  I ,  where  I  is  the  identity  matrix. 

Let  px  '  be  the  last  row  of  P  and  be  defined  by 

^  I  l©M 
I  I  PK  I  I  *  =  1 

Note  that  this  can  be  done  because  an  orthonormal  tranformation  can  always  be  completed. 
Now  let  X  -  PZ  and  calculate 

-  Pk  ^ 

- 

2j  r  : 


j=o  I 


u  ;  e 


u  10:  '  ^ 

^  ei  , 


Therefore  =  0  random  variable,  and  the  space  is  K-l  vice  K  dimensional. 
Ne.xi  calculate 

Cov(X,X)  --  Cov(PZ.PZ)=  E{PZ){PZ)' 

E(PZZ  P  ')  r  P'E{  ZZ  '  yp  ' 

P{('ot  [Z.Z))P'  -  P{1  — — -\P 


pp  ■  p 


-r  -  --  I  P  p,  p,  ■  p 


D  O'  It’  0 

P  Px  Pk  P  -  0  }  i 


Therefore, 


!• 


Cov{X.X)^  [  0  6J 


Since  X  =  PZ  implies  Z  =  P  '  X ,  and 

V  =  (Z  Z)  =  (P'Z)'(/»  X)  =  (X  PP  Z)  =  X  X, 


so  V  has  same  distribution  as  |  |  Z  |  |  ’. 
Then, 


Jf  K-l 

y«o  i-i 


Therefore,  V  is  distributed  ;^(ir-i). 


To  show 


W  =  UKV  =  .  :  T  I  1 


t/2-rV'  =  f/2+  :  I  Y 


=  vK(y-  .  ) 


!  e  I  Me:  "  .  e 


=  U^-Y  'Y~  V^-  1'* 


=  K  '  r  =  i  y 


D.  NON- CENTRAL  CHI  SQUARE 

This  section  considers  the  Non-Centr&l  Chi  Square  with  n  degrees  of  freedom.  Letting 

A 

Yi,  .  .  .  ,  Y„  be  iid  Normal  (  0,1  ),  then  X(«)  -  S  density: 


/.(«)=  — 


•.*■**■> 


0  <  n  <  00 


2*r(-|) 


But  if  the  means  are  not  all  zero  then  the  result  is  the  Non-Central  Chi  Square  random  vari¬ 
able. 

i)  With  n=l  and  letting  Y  be  distributed  Normal  (6,1)  and  V  =  Y^,  then  the  density 
of  V  is  as  follows. 


y  =  V  y 


.  /—  .  iy  \  d  I  v  * 


I  v(v)  =  \f  yW^)  ^  I  y[- 0  <  ti  <  oc 

1  1  i  -{l-vV  -bA 

=  V<  .  -(-«■'  / 

2\A>  \/2ir 

1  I  (  -i|r*I©v/7  -  eA 

=  -—7=  v«  -t-  e  *  I 

2\/v 

=  — U  -4= 

2\/  V  v'^  '•  I 

_  ^  _  Z 

=  — e  *  — ^  e  *  Co*h  (©n/iT  ) 


^  e  *  e  ^  ^  ^  f 

vu  (2")! 


1  2  2  1  <  6 


2  2 
e  e 


-4=E(^)"2" 

V  V  „.0  2  n  . 


n !  1 


-±-  t  =_Lvi - L(.^)"  ’ 

v2;r'  V  ••  „‘^o  2'  (2..)' 


E/’JA),."  ^c 

„.o  \  (2ri  )' 


Where 


2 


"  7  T 

V  ’e  * 


P»(A)  -  c 

n  : 


and  is  the  constant  for  the  central  ^  density  function. 

ii)  The  case  of  general  "n"  can  be  managed  by  transforming  all  of  the  non-centrality  into 

f| 

the  first  variable.  As  before,  V  =  2  • 

•  =■1 

Let  P  be  an  orthonormal  transformation  {P  'P  =  /)  and  having  the  first  column  pj  : 

I  I  e|  I 

Note  that  p  i  p  i  =  - -  =  1  and  that  an  orthonormal  transformation  can  always  be 

lien* 

completed. 

Then  let  T  P  '  Y 

and  hence  V=Y'Y=T'P  ’  PT  =  T'T=  \  \  T  \  \  ^ 

and  T i,  .  .  .  ,  Tn  are  independent  (i.e.,  since  T  is  a  normal  random  vector,  then  one  needs 
only  compute  the  covariance  matrix). 


Cov(T,T')=  E(TT')  -  E[T)E(T  ■) 


Now  for  i=l 


For  i  >1 


E[T,)  =  £(p,'F)  =  — j-£(K)  =  ,  =  i  i  e, 


£(7,)=  £(p/F)=  p,£(K)=  p.e:.  ©  p.p,=  0 


Hence  E(T  ')  =  ( !  '  8]  ;  ,0,  .  .  .  ,  0)  and 

E(TT')=  E{P'YY'P)=  P'E(YY')  P 

=  P  '\Cov(Y,Y')  ^  E(Y)E(Y')]  P 

=  P  1/  +  ©©  I  P 


=  P  'P  ^  P  '©©'P 


Using  P  P  =  I  and 


©  P  =  (  i  ©  .0,  .  .  .  ,0)  =  £(r 


Hence. 


Cov(T.T')  =  I  -  E{T)E{T  ')  -  E{T)E{T  I 


and  independence  follows. 

Since  E(T,)  -  0  for  i  1.  it  follows  that  the  Non-Central  \(n)(A)  where 
^  ^  is  the  convolution  of  and  a  Central  .  This  is  treated  in  the 

next  section. 


E.  DISTRIBUTION  of  W  and  EXPECTATION  of  W' 

This  section  illustrates  the  characteristics  of  the  distribution  "W".  which  is  the  convolution 
of  /yj(*  )  (non-central  chi  square)  and  the  density  of  -i)  -  //t  i(“  )  ■  noted  that  W  is  a 

non-central  chi  square  random  variable  with  K  degrees  of  freedom  and  non-centrality  parameter 

A. 


/»'(«')=  r  /»»(«-«')  /k -I  (»)  'f*' 

•h 

00  ^ 

=  r  f)  l-^)  /if-I  («’)  ‘f*' 

11*0 

00  ** 

=  S  (-')  /  /ja-ll**'-*')  fh-l  («')  ‘ft’ 


=  S  fin^K  (“') 


Desire  to  derive 


El  — 
IV 


oc  ^ 

r  -  E  f*-  (A) 

n=0 


oc 

=  E  (‘"I'fw' 


Lei  r  -  2n  -  k 


—  1  — 

Now  /,(«')  =  — ^ e  * 


2  =  r(-) 
2 


«  ^  n  —  »  ^ 


2-2 


C  ^  rfu* 


2^ny)  ” 


oc  00  r  ^  > 

E  (^)  —  f  ®  ‘  du 

n  -  C  ”  r  *  0 

2  1(7) 


Lei  11  -  —  .  (hen  du  -  —  du 

2  2 


IJ-  -  E  f’l.  -7-^ -  /  (2«)  '  “2 


2*n7)  ” 


wx»w 


F.  LOSS  AND  RISK  FUNCTIONS 


This  section  evaluates  the  Loss  and  Risk  functions.  Considering  the  column  vector  Y.  which 
is  distributed  Normal  [Q.I],  these  functions  are  as  follows. 

LOSS  =  L  (a  ,©)  =  (a  -©)(o  -©) '  =  .  a  -  © :  ^ 

Risk  is  expected  loss,  and  for  the  maximum  likelihood  estimator,  a  =  Y  and 

RISK  =  E[  L(y,©)  !  =  E(Y-  ©)'  (Y-  ©)=£(!  !  Y-  ©!  \  ^)  =  K 


Consider  another  function  of  Y.  Use  I  1  Y  1  i  *  =  Y  '  Y  =  ^  which  is  distributed 


>  =  o 


Non  Central  Y(A')-  See  section  D  for  the  properties  of  the  Non  Central  Chi  Square  Distribution. 
The  James-Stein  estimator  has  structure; 


a(Y)  =  (1- 


;  Y 


)y 


It’s  Risk  function  is  derived  as  follows: 

b 


RISK  =  E\\  ;  (1 - 

=  £ ! ,  Y  -  © 


Y  I 


-  )  J'-  ©I 


bY 


i'i  I 


1  2 


=  £  ( Y-  ©  - 
=  £|(Y-  ©) 


bV 


y  * 

bY 


( Y-  ©  - 


6Y 


)1 


(Y  -  ©)  - 


6Y 


£;(Y  -  ©)'  (^'  -  ©)  - 


I  Y  ^ 
Y'  {Y-  ©)  - 


=  £ 


y  -© 


Y  :  ’ 

-^Y'  (Y  -  ©)- 


Y  -  ©) '  V  - 


6- 


y 


=  £ 


Y-© 


-2£ 


Y  '  ©):-fc*£ 


K  -  2E 


>  '  (  Y-  ©)'^6*£' 


X(K  li-^) 


Let 


EVAL  =  £| 


Y  ©)  ■  Y 


) 


C-)  £(^) 


where  U  and  W  are  defined  in  section  B  of  this  appendix. 


U 


See  section  G  in  this  appendix  for  an  evaluation  of  E{-—)  .  I  sing  that  result  now  leads  to 

H 


the  following; 


EVAL=  10,  *f:  P,(A)  — ^ 

t  =  0  n  I J 


Where  A  =  ^ — ,  and  /*y(A)  =  e 


EVAL=  1-A  S  P,(A)  — ^ 
>*o  * 


-  ,  V  2  J+1 

>•  ^r-2y  3+1 


Let  f  =  y^l  and  note  that  the  following  expression  in  the  summation  is  zero  when  i  =  0. 
therefore  the  index  can  be  initially  zero. 


-  E/.wc-ir^i 


"  E  c,  (A|( 


K  -2i  -2-2t\ 


K  ~2i  -2 


K  -2J  2 


Where  J  i.«  distributed  Poisson  (A)  and  Prob  (J--i)  =  P,  (A). 


Therefore. 


RISK  -  K  2b  {E  -2]E  — - - -  -b‘E  — - - - 

2y*A-2  2J  ~  h  2 


Choose  b  to  minimize  Risk  .  let  5  =  £? 


2J  ^  K  -  2 


d  RISK 
d  b 


2S(h  2)  -  2bS  --  0 


b-RlSK 


2S  •  0 


.Solution  i.s  h  K  2  ihereforc 


RISK  K  2{K  2]-f:  — - - - iK  2]-E  - 

2.1  ■  K  2  i 


2J  -  A  -  2 


K  -{K  -2fE 


2J  ■  K  2 


Now  Ov  E  - ; - - -  ■  — - - 

2.1  h  2  h  2 

Hence  2  ?  RISK  K  . 


^Vu>-w*vLVw  •  V  r_k  Tj.  •ja.xVk.y.^r 


G.  EVALUATION  OF  E(  ) 


U 


The  goal  of  this  section  is  to  evaluate  E  {  )  . 


First,  find  the  joint  distribution  of  U  and  W. 
fu,v(*t,v)  =  /v{»)fv(«) 


1  -jc--  I  ieM)>  1 _  ^  ^ 


\/2v 


K-  1 


t;  *  c  * 


0  <  »  <  oo 


2  *  n-^) 


Change  of  variables; 


2  ,  2  . 
«D  =  U  +  U,  V  =  W  —  U  ,  — —  =  ] 


dw 


0<l’<OO  =>  0<tl,-U*<OO  =>  U*<ID<OC 

1 


,  ,  ,  1  -{(•-!  lei  D’ 

/  t/,W  («  W  )  =  — ^  «  y  7 

\  2jr  *  ~  ' 


K  -  i  1  ,  I, 

.2\  2  ^2'  ^ 


2  2  — 1) 

2  ' 


(u;  2  g 


“00<u<00  u  <  w  <  oo 


Since  -  -y  («  -  I  i  ©I  I  )*  =  -  y(«*-2tt  |  1  © |  |  +  |  |  © |  |  *)  it  follows  that 


f  V  ,w  (u  .w  ] 


—  K-  I 

2  *  — i-) 

2 


,  -  ■L(.-2«  M,ei  I  *l  lei 

(w-tt*)  *  e  * 


Therefore, 


-r  «  ’ 


K-  I 


00  « 

2 


J 


v/27r  2  *  r( 


A  -  1  .  ni 


/  ^(.-u 

—  U’ 


K-  I 


^1  ^  c  ® 


du 


dw 


Where  the  range  is  u 


u  ■>•:  u  V  «•  . 


Let  t  -  — .  dt  =  — ^=-  ,  and  -  1  <  f  <:  +1.  Continuing 

\  U'  \  « 


—  e 


K~  1 


X  «  I 

r  e  •  f  — L-  u 


A'  -  s  K  -  i 


(1-  t 


2\  2  ,  e  (\  u' 


—  00<  tl 


w  dt\  dw 


Let  t  =  2y+l. 


'  W 


-  v;""’  f  [1- 

2  ' 


>+»  (i_t 


Now  the  inner  integral  is  a  Beta  i  j  +-^  ,  ^  +l)  =  fleto  !;>—  ,  - ^1,  and 

2  2  2  2 


Betaij^-  ,  — 
'  2  2 


ro+i+f) 


£(  —  )  ^ 
'  W  ' 


{Me..* 


+  1  ro--!)  «  ,K_ 


«  " _ n--— 1  ^  2'  p 


The  integral  is 


I’(j'y)  yK  ^ 

- Y~  -  2  *r|y-^— ).  It  follows  that 


t  2  ' 


4-  e  ’ 


2  riJ^^ _ L 


v'2jr  2  ^  r(-^^ - -) 

2  ' 


.r{— — L)  ^2^^  *  — — — — 
2  (2y-l)! 


^  JT 
2 


^(2j-l)(2y-l)  (3)(l)r(|) 


Note  further  that 


(2y-l)’  (2j  •  1)(2;  )(2y  1)  ■(31(2)(1) 

-  (2j  -l)(2y  1)  •  (5)(3)(1)  (2,)(2)-2)  (6)(4,(2| 

-  n  (2'*«)  2^  (;!) 


it^//,V  «  V  >'/>v 


H.  INVERSE  SINE  TRANSFORM 

Method  (1): 

©,  =  \/n  ^  25in''  v'p 

Method  (2); 

©j  =  [5tn“'  (2p  -1)  +  j 

z> 

Method  (2)  was  used  because  it  was  computationally  faster  for  the  computer  to  execute 
(2p-l)  vice  .  Proof  that  both  methods  are  the  same  is  as  follows.  First  step  is  to  reduce  the 
right  hand  side  of  each  equation  to  the  elements  within  the  brackets. 

Sin  (  — ^  )  =  Sin  [25tn  ‘  y/^  ; 

V  n 

=  2\'^  Cos'Sin~*y/p  ; 

=  2v^  y/l-p 

=  2\/p  (1-p  ). 

Sin  (  — :==■  )  =  Sin  15in‘'  (2p  1)  -  — 

V  n  2 

=  (2p  -1)  Sin~*  (2p  -1) 

-  (2p-l)(0)  -  (l)v'l-{2p-l)* 

-  V  1  -(4p'-4p  -  1) 


^  2v  p  (1-p  ). 


~r  =  Sin 


©j 

— implying  0, 


©,.  Q.E.D. 


Therefore  Sin 


0, 


I.  MONTHLY  VERSES  YEARLY  RATES 

This  section  provides  a  recommended  method  to  convert  yearly  attrition  rates  to  monthly 

rates. 

Let  r  =  yearly  attrition  rates. 

Let  =  monthly  attrition  rates. 

Assume  independence  from  month  to  month.  Thus  12  consecutive  survivor  months  implies 
a  survivor  year. 

1  -  »•  =  (1  - 

(1  -  r)'^  =  1- 

I 

=  1  -  (1  -  r 

Thus,  the  above  equation  relates  the  yearly  attrition  rates  to  the  monthly  attrition  rates. 


AiSMSU  £ 

PBG6BAHS  Al£  BDMCTIOHS 

Ibis  appendix  contains  saaple  listings  of  JCL,  FOfilBAN, 
AND  API  progcans  utilized  by  the  author  throughout  this 
research  project. 

A.  0BI61NAL  SOHHABI  DATA  PBOGBAH 

The  original  sunoary  data  from  UPttDC  is  on  a  file  in  the 
Hass  Storage  at  the  HPGS  Conputer  Center  with  a  file  name  of 
"COOmS".  In  order  to  access  the  data  from  CHS  and  transfer 
a  copy  of  the  data  set  "COON IS"  from  Hass  Storage  to  the 
HVS004  disk,  submit  the  program  "HSSCOUH  JCL  A1"  to  HVS. 
See  figure  0.1. 

The  data  set  will  nearly  fill  half  a  disk  space  of  4 
cylinders,  thus  it  is  advisable  to  get  a  temporary  work 
space  while  logged  on  in  order  to  conserve  A-disk  space. 
The  exec  "GEITEHP  EXEC  Al"  can  be  used  by  typing  "GETTEHP 
8  C".  See  Figure  D.2.  A  temporary  disk  space  of  8  cylin* 
ders  will  be  useable  with  filemode  "C".  The  files  stored  on 
this  disk  space  will  be  accessible  only  while  logged  on.  To 
view  the  files  on  this  temporary  disk  space,  type  "FLISI  *  * 
C".  This  is  similar  to  FLIST  on  the  A-disk  except  it 
accesses  the  files  on  the  temporary  C-disk. 

The  "GETHVS"  system  exec  is  used  to  copy  the  summary 
data  file  to  the  temporary  C-disk  in  order  to  access  it  from 
CHS.  The  exec  will  request  the  following  identification 
information  "S2209  COGNTS". 

The  summary  data  file  has  16093  records  with  53  charac¬ 
ters  per  record.  The  number  of  records  for  each  fiscal  year 
is  listed  in  Table  XX. 
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MSSCOUN  JCL  A 

//COUNTS  JOB  (2209.5555).'MAJ  D.  D.  TUCKER-.CLASS=B 
//•MAIN  ORG=NPGVM1.2209P 

//•  ♦•*•♦•*••♦♦♦••♦*•♦•*♦♦*  note  to  users  •••**♦•♦••♦**•♦♦*♦•***♦*• 

//•  THIS  WILL  MOVE  THE  DATA  SET  COUNTS  FROM  MASS  STORAGE  TO  THE 
//*  MVS004  DISK.  FIRST  CHANGE  THE  XXXX’S  TO  YOUR  USER  ID. 

//•  CHANGE  YYYY  TO  YOUR  PROJECT  NUMBER.  THEN  SUBMIT  THE  PROGRAM 
//♦  TO  MVS.  WHEN  THE  SYSTEM  RETURNS  A  COMPLETION  MESSAGE,  USE 
//•  THE  GETMVS  SYSTEM  EXEC  TO  COPY  THE  DATA  SET  TO  YOUR  A-DISK 
//♦  IF  YOU  WISH  TO  ACCESS  IT  FROM  CMS.  BE  ADVISED  THAT  THIS 

//*  file  IS  HUGE  AND  WILL  TAKE  UP  NEARLY  4  CYLINDERS. 

//* 

//•  USE  DSN=  S2209  COUNTS  IN  GETMVS 

//•  *«**«**««***«*«**«**«««**««•**«*•**•«*  **•***«■«,***  ***•«*«**«*«* 

//  EXEC  PGM=IEBGENER 
//SYSPRINT  DD  SYSOUT=A 
//SYSIN  DO  DUMMY 

//SYSUT1  DD  DSN=MSS.F0968.COUNTS,DISP=OLD 

//SYSUT2  DD  DSN=S2209.COUNTS,VOL=SER=MVS004.UNIT=3350, 

//  SPACE=(CYL,(2,2).RLSE).DISP=(NEW.KEEP) 


Figure  D.1  asscoOH  JCL  A1. 


B.  DATA  HANIPOLAIIOM  PB06SABS 

A  technique  utilized  by  the  author  to  reduce  the  depea> 
dency  or  niaaipulatiiig  such  a  large  data  rilCr  involved 
creatinq  a  tile  of  indices  corresponding  to  pertinent  data 
in  the  original  sunoary  data  file.  A  set  of  indices  was 
created  for  the  inventories  and  losses  identified  for  each 
fiscal  year. 


1  Indices 

The  set  indices  for  tne  inventory  can  oe  created 
using  the  following  procedures. 

•  Separate  the  sumary  data  file  by  fiscal  year  data  and 
create  a  new  file  for  each  fiscal  year.  For  example, 
•'COUWT77  DATA  C”  for  fiscal  year  1977  data  on  tne 
tewporary  C-disk. 


GETTEMP  EXEC  A 


ATRACE 

•  GET  A1  CYLINDERS  OP  TEMPORARY  199  DISK  SPACE; 
■  ACCESS  IT  AS  MODE  &2:  DSK&2  AS  DISK  LABEL. 

&IF  .A1  =  .?  AGOTO  -HELP 
AIF  AN  NE  2  AGOTO  -HELP 
AIF  A2  s  A  AGOTO  -HELP 
CP  DEFINE  T33S0  AS  199  CYL  A1 
AIF  ARC  NE  0  AEXIT  -1 
ASTACK  YES 

ALABEL  =  ACONCAT  OP  DSK  A2 

AIF  ARC  NE  0  AEXIT  -2 

ASTACK  ALABEL 

FORMAT  199  A2 

CLRSCRN 

Q  DISK 

AEXIT 

-HELP 

ATYPE  ISSUE  GETTEMP  <A  CYL>  <FILEMODE> 
ATYPE  WHERE  CYL  IS  NUMBER  OF  CYLINDERS 
ATYPE  FILEMODE  -  A  IS  NEVER  AUOWED 
AEXIT 
ATRACE  ON 


Figure  0.2  6ETTEHP  EXEC  A1. 

Execute  the  NUHBEH  exec  ty  typing  WNUMBEE  <XX>*'  where 
XX  is  the  desired  fiscal  year  (such  as  77  for  fiscal 
year  1‘>77  data).  See  Figure  D.3.  This  executive  will 
load  the  Fortran  program  **INV  FOnTBAN  Al”  whici*  reads 
the  data  from  the  file  ''COUhi’77  DATA  C”.  See  Figure 
£.4. 

Ihe  output  from  "INV  FORTRAN  Al"  will  be  sent  to  a  CHS 
file  "INVXX  ARRAY  C",  where  XX  is  the  desired  fiscal 
year  (saoe  as  above). 


TABLE  XX 

BECOBOS  PEB  FISCAL  lEAB 


Fiscal  Year  Nuaber  oz  Becords 


19  77 

2203 

1978 

2231 

1979 

2337 

19  80 

2351 

1981 

2317 

1982 

2324 

1983 

2330 

FILE:  NUMBER  EXEC  A1 


4TRACE  OFF 

*  TO  EXECUTE  THIS  EXEC  TYPE  "NUMBER  <YEAR>" 

»  WHERE  YEAR  IS  ONE  OF  77-85. 

GLOBAL  TXTLIB  VALTLIB  VFORTLIB  IMSLSP  NONIMSL  CMSLIB 
GLOBAL  LOAOLIB  VFLODLIB 
FORTVS  INV  (LVL(77)) 

*  FORTVS  81  (LVL<77)) 

CP  TERM  LINESIZE  133 

8A  =  8C0NCAT  OF  COUNT  81 
8B  =  8CONCAT  OF  INV  81 

FILEOEF  08  DISK  8A  DATA  C  (RECFM  F  LRECL  53  PERM 

FILEDEF  09  DISK  8B  ARRAY  C  (RECFM  F  LRECL  9  PERM 

8TRACE  ON 

LOAD  INV  (START 

X  LOAD  81  (START 

ERASE  INV  LISTING 

ERASE  INV  TEXT 

ERASE  LOAD  MAP 

8EXIT 


Figure  D. 3  NUMBER  EXEC  A1. 
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FILE:  INV  FORTRAN  A1 


C  PROCRAM  INVENTORY 

C  *****************1 


C 

REFORMATS 

DATA  FROM  'COUNTS'  FILE  INTO  A  TWO 

#«« 

C 

DIMENTIONAL  ARRAY  OF  INVENTORY  FOR  EACH  YEAR. 

##« 

c 

THE  OUTPUT 

IS  A  SET  OF  INDICES  FOR  MATRIX 

c 

FORMAT 1 

ION 

IN  APL. 

### 

c 

c 

INDEX  SUMMARY 

### 

c 

COL 

CHAR  INFO 

#«# 

c 

1 

1 

YEAR  (1.2...,.7)-(77.78 . 83) 

### 

c 

2.3 

2 

MOS  (00. 01...., 39) 

c 

4 

1 

GRADE  (0,1...., 9) 

### 

c 

### 

5.6 

2 

LOS  (0,1,....30) 

c 

7.8.9 

3 

INVENTORY  NUMBER 

c 

#«# 

c 

USE  " 

NUMBER  <YEAR>"  EXECUTIVE 

«## 

C  ***************************************************** 

c 

c 

INTEGER  GRADE,  MOS,  YEAR,  LOS.  INDEX,  IN.  TOT.  lYR 
C 

C  THE  DO  LOOPS  CORRESPOND  TO  THE  NUMBER  OF  RECORDS 
C  FOR  YEARS  1977  TO  1983 

C  THE  FIRST  ONE  IS  THE  FULL  DATA  SET. 

C  THE  SECOND  IS  FOR  TESTING 

C 

C  DO  200  I >1,16093 

C  DO  200  1-1.50 

C  00  200  1-1,2203 

C  DO  200  1-1,2231 

C  DO  200  1-1,2337 

C  DO  200  1-1,2351 

C  DO  200  1-1,2317 

C  DO  200  1-1,2324 

DO  200  1-1,2330 

101  READ(  8, 60 ) YEAR,  MOS, GRADE,  LOS.  I  NV,  LI ,  L2,  L3.  L4,  L5,  L6.  L7,  L8 

60  FORMAT  (412.915) 

C 

C  ***  INITIALIZE  THE  INDEX  1=77, 2-78, ....  7=83 
C 

lYR  =  YEAR  -  76 
TOT  -  100000000* I YR 
C 

IF(INV.GT.O)  CALL  SUM( I N, TOT, MOS, GRADE. LOS, I NV) 

200  CONTINUE 
C 

STOP 

END 

C 

C  THIS  SUBROUTINE  CREATES  THE  INDEX  ARRAY 
C  FOR  AN  INVENTORY  DATA  ELEMENT 

C 

SUBROUTINE  SUM( INDEX, l,J,K.L,NUM) 

INTEGER  INDEX, I.J.K.L.NUM 

INDEX=  I +(J *1000000 )+(K*100000)+(L*1000)+NUM 
WRITE  (9,500)  INDEX 
500  FORMAT  (19) 

RETURN 

END 


Figure  0.4  IN?  FOBTBiN  A1 


•  Ite  CHS  file  "INVXX  ABfiAI  can  be  read  into  an  API 

workspace  and  assigned  to  a  naaeric  variable  by 
execnting  the  user  friendly  "CMSAPL  EXEC  A1".  See 
Figures  0.5  and  0.6.  Note  that  "CMSAPL"  has  a  provi¬ 
sion  to  add  a  variable  to  an  existing  APL  workspace  or 
to  create  a  new  workspace  each  time  "CMSAPL"  is 
executed. 

The  file  of  indices  created  in  this  manner  can  be 
utilized  to  generate  matrices  in  the  APL  workspace.  The  two 
dimensional  array  for  each  fiscal  year's  inventory  is  char¬ 
acterized  by  a  length  subject  to  the  number  of  records  in 
the  summary  data  file  for  a  particular  fiscal  year  in  which 
there  is  a  nonzero  inventory  count.  The  columns  of  the 
records  in  the  index  file  is  described  in  Table  XXI 

2 .  Loss  Indices 

The  set  indices  for  the  losses  can  be  created  in  a 
manner  similar  to  the  indices  for  the  inventory  by  using  the 
following  procedures. 

•  Separate  the  summary  data  file  by  fiscal  year  data  and 
create  a  new  file  for  each  fiscal  year.  For  example, 
"L0SS77  DATA  C"  for  fiscal  year  1977  data  on  the  tempo¬ 
rary  C-disk. 

•  Execute  the  "LOSS  EXEC  Al"  by  typing  "LOSS  <XX>"  where 
XX  is  the  desired  fiscal  year  (such  as  77  for  fiscal 
year  1977  data) .  See  Figure  D.7.  This  executive  will 
load  the  Fortran  program  "LOSSES  FOfilBAN  Al"  which 
reads  the  data  from  the  file  "C0UN'i77  DATA  C".  See 
Figure  D.8. 

•  The  output  from  "LOSSES  FOBTBAN  Al"  will  be  sent  to  a 
CMS  file  "LOSS XX  ABBAX  C",  where  XX  is  the  desired 
fiscal  year  (same  as  above). 

•  The  CMS  file  "LCSSXX  ABBAX  C",  can  be  read  into  an  APL 
workspace  using  similar  methods  as  the  inventory  file 
described  above. 
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CUSAPL  EXEC  A 
ATRACE 

ATYPE  ENTER  FILENAME.  PILETYPE.  MODE  OF  PILE  TO  BE  READ  INTO  APL 
AREAD  ARCS 

AIF  AN  NE  3  AGOTO  -TELL 

STATE  A1  A2  AS 

AIF  ARC  NB  0  AGOTO  -ERROR 

ATYPE  ENTER  THE  NAME  OF  THE  APL  VARIABLE  THAT  WILL  STORE  THE  DATA 
AHEAD  VAR  AA 

ATYPE  ENTER  NAME  OP  APL  WORKSPACE  YOU  DESIRE  THE  VARIABLE  AA  IN 
AREAD  VAR  AWKS 

ATYPE  7?r?  WHAT  TYPE  OF  VARIABLE  IS  DESIRED? 

ATYPE  ENTER  N  IP  YOU  DESIRE  A  NUMERIC  VARIABLE. 

ATYPE  ENTER  C  IP  YOU  DESIRE  A  CHARACTER  VARIABLE. 

ATYPE  ENTER  Q  IF  YOU  DESIRE  TO  QUIT 
AREAD  VAR  AASK 
AIF  AASK  EQ  Q  AEXIT 

ATYPE  YOUR  INFORMATION  WILL  BE  STORED  IN  APL  WORKSPACE  AWKS 

AIF  AASK  EQ  N  AGOTO  -NUMS 

CP  TERMINAL  APL  ON 

ASTACK  )LOAD  980  CMSIO 

ASTACK  AA  'CMSREAD 

ASTACK  At 

ASTACK  A2  A3 

ASTACK  C 

ASTACK  )WSID  MYSPACE 
ASTACK  )SAVE 
ASTACK  )OFP  HOLD 
EXEC  APL 
AGOTO  -SKIP 

-NUMS  CP  TERMINAL  APL  ON 
ASTACK  )LOAO  990  CMSIO 
ASTACK  AA  'CMSREAD 
ASTACK  A1 
ASTACK  A2  A3 


Figure  D.5  CHSiPL  BXEC  A1. 
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&STACK  N 

4STACK  )WSID  MYSPACE 
iSTACK  )SAVE 
itSTACK  )OFP  HOLD 
EXEC  APL 

-SKIP  icTYPE  ■>'>??  IS  THIS  A  NEW  WORKSPACE?  Y  OR  N 

jcREAD  VAR  &NEW 

JcIF  ANEW  EQ  Y  &GOTO  -JUMP 

CP  TERMINAL  APL  ON 

JcSTACK  )LOAD  JcWKS 

ASTACK  )PCOPY  MYSPACE  *A 

iSTACK  )SAVE 

iSTACK  )DROP  MYSPACE 

4STACK  )OFP  HOLD 

EXEC  APL 

JcEXlT  98 

-JUMP  CP  TERMINAL  APL  ON 
iSTACK  )W3ID  iWKS 
JtSTACK  )PCOPY  MYSPACE  JeA 
JcSTACK  )SAVE 
ASTACK  )DROP  MYSPACE 
iSTACK  lUFP  HOLD 
EXEC  APL 
AEXIT  99 

-TELL  ATYPE  YOU  HAVE  ENTERED  TOO  MANY  OR  NOT  ENOUGH  ENTRIES 
ATYPE  FOR  THE  FILE  TlUT  YOU  WANT  TO  BE  READ  INTO  APL. 

ATYPE  YOU  NEED  TO  BEGIN  AGAIN 
ATYPE 

ATYPE  ENTER:  CMSAPL  ' 

AEXIT  too 

-ERROR  ATYPE  A1  A2  A3  DOES  NOT  EXIS  ON  YOUR  JcCONCAT  OF  A3  -DISK 
ATYPE  CHECK  YOUR  FLIST  AND  THEN  BEGIN  AGAIN 
ATYPE  BY  ENTERING  CMSAPL 

AEXIT  tOI 


Figure  D.6  caSAPL  EIEC  A1  (COMT.) 


^  m 

y 

V. 


0. 


TABLE  XXI 

IHVEIITOBX 

IMDEX  SOHtlABY 

Column 

Characters 

I  nformation 

1 

1 

YEAR  (1=1977... 

2 

MOS  (Jo...  39) 

1 

GRADE  (0...9) 

7?&% 

2 

LOS  (0 _ 30) 

3 

INVENIORY  NU:iB£H 

FILE:  LOSS  EXEC  A1 


STRACE  OFF 

*  TO  EXECUTE  THIS  EXEC  TYPE  "LOSS  <YEAR>" 

X  WHERE  YEAR  IS  ONE  OF  77-85. 

GLOBAL  TXTLIB  VALTLIB  VFORTLIB  IMSLSP  NONIMSL  CMSLIB 
GLOBAL  LOADLIB  VFLODLIB 
FORTVS  LOSSES  <LVL(77)) 

X  FORTVS  81  (LVL(77)) 

CP  TERM  LINESIZE  133 
8A  =  SCONCAT  OF  COUNT  81 
8B  =  8C0NCAT  OF  LOSS  81 

FILEDEF  08  DISK  8A  DATA  C  (RECFM  F  LRECL  53  PERM 
FILEDEF  09  DISK  8B  ARRAY  C  (RECFM  F  LRECL  9  PERM 
8TRACE  ON 

LOAD  LOSSES  (START 
X  LOAD  81  (START 
ERASE  LOSSES  LISTING 
ERASE  LOSSES  TEXT 
ERASE  LOAD  MAP 
8EXIT 


Figure  D.7  LOSS  SXEC  A1. 
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PROGRAM  LOSSES 

***  REFORMATS  DATA  FROM  'COUNTS'  FILE  INTO  A  TWO  *** 
***  DIMENTIONAL  ARRAY  OF  LOSSES  FOR  EACH  YEAR.  *** 
***  THE  OUTPUT  IS  A  SET  OF  INDICES  FOR  MATRIX  *** 

***  FORMATION  IN  .APL.  *** 

***  *** 

***  INDEX  SUMMARY  *** 

***  COL  CHAR  INFO  •** 

***  1  1  YEAR  ( 1,2. . . . .7)  =  (77, 78 . 83)  **• 

***  2,3  2  MOS  (00,01, ... ,39)  *** 

***  4  1  GRADE  (6,1,.... 9)  *** 

***  5.6  2  LOS  (0,1, .. .,36)  *** 

***  8,6,10  3  NUMBER  0^•  LOSSES  *** 

***  *** 

***  USE  ”LOSS  <YEAR>**  EXECUTIVE  *** 

***************************************************** 

INTEGER  GRADE, MOS,  YEAR,  LOSS,  LOS,  lYR 


INDEX  SUMMARY 
INFO 


THE  DO  LOOPS  CORRESPOND  TO  THE  NUMBER  OF  RECORDS  FOR  YEARS  77-83 
THE  FIRST  ONE  IS  THE  FULL  DATA  SET,  THE  SECOND  IS  FOR  TESTING 


DO  200  1= 
DO  200  1= 
DO  200  1= 
DO  200  1= 
DO  200  1= 
DO  200  1= 
DO  200  1= 
DO  200  1  = 
DO  200  1= 


1,16093 

1-50 

1,2203 

1,2231 

1,2337 

1,2351 

1,2317 

1,2324 

1,2330 


**  LI  -  L8  IS  LOSS  TYPE  CODE  1 


READ (8,60) YEAR . MOS , GRADE , LOS , INV, LI , L2 , L3 , L4 , L5 , L6 , L7 , L8 
format  (412, 9IS) 


***  INITIALIZE  THE  INDEX  1=77,2=78 . 7=83 


lYR  =  YEAR  -  76 
LOSS  =  100000000* lYR 


IF(L1.GT.0 
IF  L2.GT.0 
IF(L3.GT.O 
IF  L4.GT.0 
IF  L5.GT.0 
IF  L6.GT.0 
IF^  L7.GT.0 
IF  L8.GT.0 
CONTINUE 
STOP 
END 


CALL  SUM(IN,1 
CALL  SUM  IN,  I 
CALL  SUMi  IN,1 
CALL  SUM  IN,  I 
CALL  SUM  IN,  I 
CALL  SUM  IN,  I 
CALL  SUMi  IN,  I 
CALL  SUM  IN,  1 


,  MOS, GRADE, 
,  MOS , GRADE , 
.MOS, GRADE, 
,  MOS , GRADE , 
,  MOS, GRADE, 
,  MOS , GRADE , 
,  MOS, GRADE, 
,  MOS, GRADE, 


LOS, 1, LI) 
LOS,2,L2' 
LOS , 3 , L3 
LOS, 4, L4 
LOS,5,L5 
LOS, 6, L6 
LOS , 7 , L7 ' 
L0S,8,L8 


THIS  SUBROUTINE  CREATES  THE  INDEX  ARRAY  FOR  LOSS  DATA  ELEMENT 


SUBROUTINE  SUM( INDEX, I , J, K, L, M, NUM) 

INTEGER  INDEX, I, J.K.L.ft.NUM 

INDEX=  I ♦ ( J* 1 OOOOOOO W ( K* 1000000 ) ♦ ( L* 10000 ) + ( M* 1000 ) ♦NUM 
WRITE  (9,500)  INDEX 
FORMAT  (110) 

RETURN 

END 


Figure  0.8  LOSSES  FOBIBAN  A1. 


Ihe  file  of  indices  for  ttie  losses  is  similar  to  the 


inventory  file.  The  two  dioiensional  array  ror  each  fiscal 
year*£  loss  is  character! zea  by  a  length  subject  to  the 
nuuber  cf  records  in  the  suiimary  data  tile  for  a  particular 
fxscal  year  in  whicu  there  is  a  nonzero  loss  count  for  each 
loss  type.  It  should  be  noted  that  a  record  in  the  original 


suanary  data  aay  have  a  naaerical  count  for  aore  than  one 
type  cf  loss.  The  coluans  of  the  records  in  the  index  file 
is  described  in  Table  XXII. 


TABLE  XXII 
LOSS  mOEX  SOHIIARI 


Coluan 

Characters 

Inf or nation 

1 

1 

YEAR  (1=1977.. .7 

2 

EOS  ((70...  39) 

1 

GRADE  (0...9) 

2 

LOS  (0...30) 

1 

LOSS  TYPE  (1... 8) 

8,9,10 

3 

NUhBER  CF  LOSSES 

C.  fONCTIOBS 

Nuaerous  APL  functions  were  utilized  in  this  project  for 
data  lacipulation  and  execution  of  calculations  pertaining 
to  the  processes  under  evaluation.  The  following  is  a 
discussion  of  the  aost  utilized  functions  and  their 
purposes. 

1 .  Inventory  Functions 

The  workspace  ''IN¥TOT"  was  utilized  to  retain  the 
inventory  index  aat rices.  The  variables  assigned  to  each 
fiscal  year's  inventory  index  aatrix  was  labeled  "INVXX", 
wnere  "XX"  is  the  applicable  fiscal  year  (e.g.,  77  for 
fiscal  year  1977)  . 


a.  Function  GBTINV 

The  function  to  create  the  inventory  iLatcix  for 
each  desired  year  is  "GETINV*.  See  figure  D.  9.  "GETINV" 
utilizes  the  function  '*lNVhATX  IMVXX",  which  interprets  the 
inventory  index  natrix  "INVXX"  and  creates  the  inventory 
matrix  for  fiscal  year  "XX”.  See  Figure  0.10.  The  resul- 
tant  matrix  is  "IXX"  for  fiscal  year  "XX".  The  function 
"INVHAIX"  could  create  a  matrix  of  the  following  dimension 
7x40x10xJ1  tor  7  years,  40  {!OS*s,  10  grades,  and  31  LOS's. 
However,  due  to  limited  workspace,  the  dimension  of  40x31x10 
for  40  MOS*s,  31  LOS*s,  and  10  grades  was  commonly  utilized. 


~ 

'  .  . 

V  GETINV 

T 

«  GET  THE  NUMBER  OF  YEARS  OF  INVENTORY  DESIRED. 

'2' 

a  COMMENT  OUT  YEARS  NOT  USED  IMMEDIATELY. 

;3 

I77«-1NVMATX  INV77 

"4 

178-INVMATX  1NV78 

‘5 

I79«-INVMATX  INV79 

’6 

leO^INVMATX  INV80 

'7 

ISI^INVMATX  INV81 

'8 

a  182-INVMATX  1NV82 

‘9 

a  I83«-INVMATX  INV83 

:io] 

■  SHAPE  OF  177  IS  • 

[11] 

□  -pl77 

Figure  0.9  APL  Function  GETIN?. 


b.  Group  Inventory  Functions 

There  is  a  separate  function  to  make  a  combined 
matrix  for  each  of  the  NOS  groups.  Listed  in  Table  XXIII 
are  the  groups,  the  function  used  to  create  the  desired 
matrix,  and  the  shape  of  the  resultant  matrix.  Figure  D.11 
is  a  representation  of  the  programing  technique  used  ror  the 
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— 

7  Z-INVMATX  X:A;B;C;D;E;F;I:J 

'1 

••■MATRIX  INDEXING  (YEAR.MOS.GRADE.LOS.INVENTORY)  ' 

'2 

Z-(40  31  10)p0 

3 

«  X  IS  AN  ARRAY  OF  INDICES  FOR  A  YEAR  OF  INVENTORY  DATA. 

'4 

«  INDEX  =  9  CHARACTERS  PER  ARRAY  ELEMENT 

■5 

«  INDEX  NUMBER  CHARACTERS  DESCRIPTION 

'6 

n  1  1  YEAR  (1...7)  FOR  (77.. .83) 

7 

«  2.3  2  MOS  (00,01.. ...39) 

8 

A  4  1  GRADE  (0.1 . 9) 

■9 

A  5.6  2  LOS  (00.01 . 30) 

’10 

A  7,8.9  3  INVENTORY 

‘1 1 

a 

’12 

l-pX 

13 

J«-l^9 

'14 

LOOP:-‘(J=0)/OUT 

'15 

A  THIS  INDEX  IS  NOT  USED  BECAUSE  WITH  THIS  ADDED 

‘16 

A  DIMENSION.  YOU'LL  GET  A  FULL  WORKSPACE. 

‘17 

A-»(1fX) 

'18 

A 

'19 

B-1+(A(2tX-(UX))) 

’20 

C«-1+(*(1tX-^(2AX))) 

’21 

D^H-^(2tX^0AX))) 

'22 

E-A(3rX-(2AX)) 

‘23 

2[B;D;C]-E 

'24 

X-(3lX) 

25 

J-J-1 

'26 

-»L00P 

27 

0UT:'FINISHED  —  SHAPE  OF  MATRIX  IS  ' 

[28] 

PZ 

i _ 

. .  . . .  . . -  -  .  ■■■_■  _  ■_ 

Figure  D.  10  iPL  Fanction  IHVHATX. 

Ground  Combat  NOS.  Ihe  technique  for  the  other  groups  was 
similar. 

These  functions  call  the  function  "GETMuS  fl", 
which  creates  the  central  data  inventory  for  the  desired 
fiscal  year  for  a  particular  NOS  of  "YY".  In  order  to  get  a 
comoat  Service  Support  NOS  inventory  matrix,  the  “GEINOS  YY" 
function  can  be  utilized  separately  to  get  the  NOS  "yy" 
matrix.  The  function  *'GETN0S'*  uses  the  glonal  variables  of 
"liY"  and  “LXX*',  for  tne  inventory  matrix  and  loss  matrix 
respectively,  for  fiscal  year  ”XX".  See  Figure  D.  1^. 
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7  GETGC;K 

;il  «  THIS  GETS  THE  INVENTORY  MATRIX.  FOR  THE  GROUND 
'2]  «  COMBAT  OCCUPATION  GROUP. 

;3]  B  SHAPE  IS  (MOS;YEARS;LOS;GRADE) 

4]  n  •*  CHECK  THE  2NO  INDEX  FOR  THE  CORRECT 

5]  »  *•  NUMBER  OF  YEARS  DESIRED. 

;6]  GC^(3  4  31  10)p0 

■7]  a  THIS  GETS  THE  MOS;  INFANTRY 
’8]  GC[1;;;]^GETM0S  3 

■9]  a  THIS  GETS  THE  MOS;  ARTILLERY 
;iO]  GC[2:::]«-GETMOS  5 

’1 1]  a  THIS  GETS  THE  MOS;  TANKS  AND  AMPHIB 
T2]  GC[3;;;]-GETMOS  10 


Figure  0.11  APL  Function  GET6C. 


7  Z^GETMOS  X;J;Y 

T]  a  GET  THE  CENTRAL  DATA  NUMBER  OF  INVENTORY 
[2]  a  FOR  THE  YEARS  DESIRED  FOR  A  PARTICULAR  MOS. 
■3]  a  X=  MOS  DESIRED.  CHANGE  THE  INDEX  FOR  Z  IF 
■4]  a  YEARS  ARE  NOT  SEQUENTIAL. 

■5]  J^X-l-1 

“ej  Z«-(4  31  10)p0 

y]  Z[1  ;;]^((Y-(l77[j;.]+l78[J;;])^2)r(L77[j;;])) 

8]  2[2;;]-((Y-(l78[J;;]+l79[J;;])^2)f(L78[J;;])) 

>1  Z[3;.]-((Y-(t79[j;;]+180[J;.])T2)f(L79[j;;])) 

;iO]  Z[4;;  -((Y^(l80[j;;  +I81[J;;  )^2)r(L80[J;;])) 

;il]  aZ[5;;  »-((Y^(l8l[j;;  +l82[J;;j)^2)f(L8 1  [J;;])) 

;i2]  «Z[6;;  -((Y-(l82[j;;  +I83[J;;])^2)[(L82[J;.])) 


Figure  0.12  APL  Function  GETHOS. 


To  create  the  in?entory  matrix  used  ror  the 
James-Stein  tunction  the  lollowing  functions  were  used: 

•  GEING  -  creates  the  ”NG**  matrix  tor  Ground  Comtat.  See 
Figure  D.13  for  an  example  of  the  program  utilized  to 
create  the  the  Ground  combat  inventory  matrix.  The 


TiBL£  2XII1 
IITEHTOBT  FOMCTIOHS 


Occapatioa  Group  Puaction  Mane 


Aviators 
Coabat  Support 
Ground  CoBoat 


GEIAV 

GETCS 

GETGC 


Sha  pe 

iYears,i<us ,  grade) 

OS, Years , LOS, grade) 
MOS, Years ,LOS, grade) 


COBcat  Support  inventory  matrix  was  created 
sioilar  manner. 

•  GElhC  -  creates  the  ••NC  aatrix  for  Coabat  Support 


V  GETNG 

;i]  «  THIS  CREATES  THE  INVENTORY  FOR  GROUND 
2]  <*  COMBAT.  YEARS.  AND  GRADE  ELEMENTS. 

■3]  <*  NG(X).  X=0-8  FOR  WO  TO  LT  COL. 

■4]  «  NG0^(3  ^  2)>s(GC[;;;^3 
■5]  •  NGl-(3  1  2)5!(GC[;;;2] 

‘6]  fl  NG2^(3  1  2)^GC[:;;3j 
■7]  «  NG3»-(3  1  2)^GC[;;;4] 

8]  «  NG4-(3  1  2)>^GC[;;;5j 
■9]  NG5«-(3  1  2)^siGC[;;;6] 

’10]  «  NG6*-(3  1  2)>HGC[;;;7] 

'll]  fl  NG7^(3  1  2)5s(GCC;::8] 

’12]  NG3-(3  1  2)!stOC[;,;9] 


Figure  0.13  APL  Function  GETNG. 


2.  Logs  functions 

lae  workspace  "ilATBiX"  was  utilized  to  retain  the 
loss  index  matrices.  The  varianles  assigned  to  each  fiscal 
year's  less  index  matrix  was  labeled  "LOSSXX",  where  "XX"  is 
the  applicarle  fiscal  year  (e.  g.,  77  ror  tiscal  year  1977). 


a.  Function  GETLOSS 

The  function  to  create  the  loss  aatrix  for  each 
desired  year  is  "GETLOSS**.  See  Piyure  D.  14.  "GETLOSS" 
utilrzes  the  function  "N&TBIE  LOSSXX'*,  vhicb  interprets  the 
loss  index  natrix  "LCSSXX"  and  creates  the  loss  matrix  for 
fxscal  year  "XX".  The  resultant  natrix  is  "LXX"  for  fiscal 
year  '*XX**.  The  function  "flAIfilX**  could  create  a  matrix  of 
the  following  dimension  7x40x10x31x8  for  7  years,  40  NOS*s, 
10  grades,  31  LOS*s,  and  8  loss  types.  However,  due  to 
limited  workspace,  the  dimension  of  40x31x10  for  40  NOS's, 
31  LCS's,  and  10  grades  was  commonly  utilized  when  the 
losses  were  aggregated.  Each  loss  type  could  be  examined 
witn  the  same  dimensicced  natrix.  See  Figure  0.15. 


‘1 

7  GETLOSS 

«  THIS  GETS  THE  L(XX)  MATRIX  FOR  YEARS 

’2 

A  XX=77 . 83. 

'3 

L77«-MATRIX  L0SS77 

■4 

L78-MATRIX  L0SS78 

■5 

L79-MATRIX  LOSS79 

‘6 

L80^  MATRIX  LOSS80 

■7 

aL81 -MATRIX  L0SS81 

8 

«L82-MATRIX  L0SS82 

[9] 

-L83-MATRIX  L0SS83 

Figure  0.14  A  PL  Function  GETLOSS. 


b.  Group  Loss  Functions 

There  is  a  separate  function  to  make  a  combined 
matrix  for  each  of  tne  NOS  groups.  Listed  in  Table  XXIV  are 
the  groups,  the  function  used  to  create  the  desired  matrix, 
and  the  shape  of  tne  resultant  matrix.  See  Figure  D.16  for 


7MATRIX[D]7 

7  Z'-MATRIX  X:A:B:C:D;E:F:l:J 

■MATRIX  INDEXING  (YEAR.MOS.GRADE.LOS.LOSS  TYPE.NUM  OF  LOSSES)  ‘ 
USE  FOLLOWING  DIMENSION  FOR  Z  IF  SHAPE(MOS,GRADE.LOS.LOSSTYPE) 
Z-(40  10  31  8)/j0 

USE  THE  FOLLOWING  DIMENSION  FOR  Z  IF  AGGREGATING  THE  LOSSES 
Z-(40  31  10)p0 

USE  THE  FOLLOWING  DIMENSION  FOR  Z  IF  SHAPE  (YEAR.LOSSTYPE.GRADE) 
Z-(l.8,10)p0 

USE  THE  FOLLOWING  DIMENSION  FOR  Z  IF  SHAPE  (YEAR.LOSSTYPE.MOS) 
Z*-(1.8.40)p0 

I  USE  THE  FOLLOWING  DIMENSION  FOR  Z  IF  SHAPE  (YEAR.LOSSTYPE.LOS) 
Z-(1.8.31)p0 

I  X  IS  AN  ARRAY  OF  INDICES  FOR  A  YEAR  OF  LOSS  DATA. 

I  INDEX  =  10  CHARACTERS  PER  ARRAY  ELEMENT 


INDEX 

NUMBER  characters  DESCRIPTION 

1 

1 

YEAR  (1.2 . 7)  FOR  (77.' 

2.3 

2 

MOS  (00.01 . 39) 

4 

1 

GRADE  (0.1,. ...9) 

5.6 

2 

LOS  (00,01 . 30) 

7 

1 

LOSS  TYPE  (1.2 . 8) 

8.9,10 

3 

NUMBER  OF  LOSSES 

l*-pX 

LOOP:-‘(J=0)/OUT 

«  THIS  INDEX  IS  NOT  USED  BECAUSE  WITH  THIS  ADDED  DIMENSION  WS  FULL 


[39' 

[40 

[41' 

[4-2 

[43 

[44 

[45 

[46 

[47' 

[48' 


A*-*(1tX) 

a 

B«-1+(»(2tX^(UX))) 

C*-1+(i(1tX^(2iX))) 

D*-1+(i(2rX*-(UX))) 

E-i(1tX.-(2iX)) 

F^«(3tX-(UX)) 

■  USE  THE  FOLLOWING  IF  YOU  WANT  SEPARATE  LOSSES 
a  need  to  CHANGE  THE  DIMENSION  OF  Z  ABOVE  TO  Z[40  10  31  8] 
a  Z[B;C.O:E]-F 

a  USE  THE  FOLLOWING  IF  YOU  WANT  THE  LOSSES  AGGREGATED 
a  Z[B;D;C]-Z[8.D;C]+F 

a  USE  THE  FOLLOWING  IF  YOU  WANT  SHAPE  (YEAR.LOSSTYPE.GRADE) 
0  Z[1;E.C]*-F 

a  USE  THE  FOLLOWING  IF  YOU  WANT  SHAPE  (YEAR.LOSSTYPE.MOS) 
a  Z[1;E;B]-F 

a  USE  THE  FOLLOWING  IF  YOU  WANT  SHAPE  (YEAR.  LOSSTYPE.LOS) 
Z[1;E;D]-F 
X«-(3iX) 

J-J-1 

-LOOP 

OUT;'FIMISHED  —  SHAPE  OF  MATRIX  IS  ’ 


Figure  0.15  APL  Function  BiTBIX 


the  Ground  Coabat  function. 


The  other  groups'  functions 


7  GCLOSS 

’1]  •>  THIS  GETS  THE  LOSSES  FOR  GROUND  COMBAT 
’2]  «  FOR  6  YEARS  77  -  82. 

'3]  »  SHAPE  IS  (MOS.YEAR.LOS.GRADE). 

‘4]  «  CHANGE  YEAR  INDEX  TO  DESIRED  NUMBER 
■5]  GCL^(3  4  31  I0)p0 
’6]  GCL[l;1:;]^L77r4;;] 

‘7]  GCL[1:2:;]-L78C4::3 

8]  GCL[1  ;3;;KL79[4;;] 

■9]  GCL[l;4:;]^L80[4;:i 

‘10]  «  GCL[1;5:;]^L81[4;;] 
n]  «  GCL[1;6;;]*-L82[4;;] 

;i2]  GCL[2;1;;]-L77j;6;;; 

13]  GGL[2;2;;]-^L78[6:; 

‘l4]  GCL[2:3::]^L79[6:: 

’15]  GCL[2;4;:]«-L80C6;: 

'16]  «  GCL[2;5;;]^L81[6;;] 

’17]  «  GCL[2;6;;]-L82[6;:] 

18]  GCL[3;1;;]^L77[11;;] 

‘19]  GCL[3;2;;]^L7a[l1;;] 

'20  GCL[3;3;:]^L79Ci1::3 

'21]  GCL[3;4;;]^L80[11;;] 

22]  •  GCL[3;5;;]^L81[11;;] 

’23]  B  GCL[3:6;:]^L82[11;;] 


Figure  D.16  APL  Function  GCLOSS. 

These  functions  use  the  glohal  variables  "LX4"« 
for  each  "LX"  fiscal  year.  In  order  to  get  a  Coaoat  Service 
Support  MOS  loss  natrix,  the  "LXX"  natrix  was  used  directly 
to  renove  the  desired  coabination  of  cells. 

To  create  the  loss  oatrix  used  for  the 
Jaues-Stein  function  the  following  functions  were  used: 

•  GElYC  -  creates  the  "YC"  aatrix  for  Combat  Support. 

•  GCIYG  -  creates  the  "YG"  matrix  for  Ground  Combat.  See 
figure  D.  17  for  an  example. 


c-r. 


TiBLE  EXIT 
LOSS  FO ECHOES 


1 


Occupation  Group 


3 .  Janes-Stein  f unctions 

The  following  notation  is  required  for  the  following 
f unct ions. 

•  Let  "i”  stand  for  LOS,  then  i=0,...,30. 

•  let  "j”  stand  for  rtOS,  (then  values  of  j's  depend  on 
which  nos  group  is  being  analized)  . 

•  let  Dij  =  Incidence  Matrix  of  nonstruct ural  zeroes, 
which  is  the  saue  for  ail  years.  Dij  =1,  if  cell  is 
xeinLer  of  the  feasible  region;  Dij  =  0,  otherwise. 


•  let  K  =  nuBlier  of  feasible  cells,  i.e. ,  sum  of  all  Dij. 

•  let  li  j  s  nuBber  of  leavers  in  cell  (i,  j) . 

•  Let  t  -  vhere  I  =  nuaber  of  years  of  data  used 

to  create  the  estiaator. 

•  let  Nij  =  Inventory  in  cell  (i,j)  =  Max 

((N(t)  ♦N(t*l))/2,I(tn  - 

Ihe  folloving  APL  functions  vere  utilized  in  the 
Jaaes-Stein  aodel.  See  Figure  D. 18  for  the  APL  listings  of 
the  following  functions. 

a.  Function  ABCSIM 

Let  ISM  =  Mij  ABCSIM  PMHij.  The  function 
"ABCSIM"  returns  the  inverse  sine  transf oraa tion  for  use 
when  the  success  probability  is  estimated  directly  (e.g.,  by 
MIMMAl)  . 

b.  Function  EIMCONV 

Let  Pij  =  Jij  BINCONV  Mi j  (t) .  The  function 

"BINCCMV"  inverts  the  inverse  sine  transf  oraa  ti  on. 

c.  Function  BINPBEP 

Let  isij(t)  *  Yij<t)  BINPBEP  Nij(t).  The  func¬ 
tion  "BIMPBEP"  prepares  the  Freeaan-Tuhey  version  of  the 
inverse  sine  transforaation  for  binoaial  data. 

d.  Function  JAMES 

Let  Jij  =  Dij  JAMES  ISi j  (t) .  Ihe  function 
"JAMES"  returns  a  Jaaes-Stein  estiaator  for  the  aeans  of  the 
aeans  of  the  cells  in  the  last  two  dimensions  of  the  input 
aatrix  while  being  screened  by  the  incidence  matrix  of 
non structural  zeroes. 

e.  Function  HIMMAX 

Let  PMHij  =  Yi j  (t)  HIMMAX  Mij  (t)  .  Tae  function 


7  X-N  ARCSIN  P 

«  X  IS  THE  INVERSE  SINE  TRANSFORMATION  FOR  USE  WHEN  THE 
«  SUCCESS  PROB.  P  IS  ESTIMATED  DIRECTLY  (E  G.  BY  MINIMAX) 
«  Y  IS  LOSSES;  N  IS  INVENTORY 
X-((0.5+N)«0.5)x-lo- 1  +2xP 

7  B^T  BINCONV  N;VO:V1 
«  INVERTS  ARC  SIN  TRANSFORMATION 

«  T  HAS  RANK  2;  N  HAS  RANK  3.  THE  1st  DIMENSION  IS  TIME. 
B<-0.5x1  +  1oT-r(0.5+(N*-+/N)-i-1rpN)*0.5 
V0«-T<(-o^2)x(N+0.5)*0.5 
VI  ^T>(o42)x(N+0.5)*0.5 
B*-V1+BxVOaV1 

7BINPREP[Q]7 
7  X^Y  BINPREP  N;P 

«  PREPS  THE  FREEMAN-TUKEY  VERSION  OF  THE 
«  ARC  SIN  TRANS  FOR  BINOMIAL  DATA 
«  Y  IS  LOSSES:  N  IS  INVENTORY 
X«— 1o-1+2xY^N+1 

X<^0.5x((0.5+N)*0.5)xX+-1o-1+2x(y+1)-^N  +  1 


7  P^O  JAMES  Z;K;N:ZB;ZBB:S;M:C 
«  P  IS  A  JAMES- STEIN  ESTIMATOR  FOR  THE  MEANS  OF  THE 
«  CELLS  IN  THE  LAST  TWO  DIMENSIONS  OF  Z  (SCREENED 
«  BY  THE  INCIDENCE  MATRIX  D  OF  NON  STRUCTURAL  ZEROS). 
S^-D  SUMSO  Z 
ZBB-(pZB)pZBB 

C-1  -(NxK-3)-f(2-K-KxN)x-r/S 
P^OxZBB+CxZB-ZBB 

7  X«-D  SUMSO  Z;SSE 

«  X  IS  SSB.SSE  FOR  THE  CELLS  OF  Z  (WHICH  HAS  RANK  3) 
fl  D  (OF  RANK  2)  IS  THE  INCIDENCE  MATRIX  FOR  THE  LAST  TWO 
«  DIMENSIONS  OF  Z.  K  IS  THE  NUMBER  OF  TREATMENTS. 
K^+/.0 

Z8*-0x(+/Z)-rN«-UpZ 

ZB8*-(+/.ZB)-rK 

X-+/.(Z-(pZ)pZBB)*2 

X-(X-SSE).SSEi-+/.(Z-(pZ)pZB)*2 


Figure  0.18  API  functions  in  Janes-stein  Algoritba. 


'*niNMAX"  returns  the  RINMAX  estioates  for  the  biuoaial.  The 


function  was  utilized  in  the  validation  procedures  to 
compare  the  James-Stein  estimator  to  the  MINnAX  estimator. 
See  Figure  D. 19  for  the  API  listing  of  this  function. 


r 

V  P*-Y  MINMAX  N;YY 

[1]  «  RETURNS  MINMAX  EST'S  FOR  BINOMIAL 

[2]  «  INCLUDES  REDUCTION  OVER  THE  1st  RANK  (TIME)  OF  THE 

[3]  «  INPUTS.  NN  (GLOBAL)  IS  CREATED  FOR  USE  IN  ARC  SINE. 

[4]  YY«-+/Y 

[5]  NN*-+/N 

[6]  P-(0.5+(0iiSYY)xYY^NN*0.5)-r1  +NN.0.5 


Figure  0.19  API  Function  NlNHAX. 

^  •  Validation  Functions 

The  following  APL  functions  are  an  example  of  the 
functions  utilized  in  the  valrdation  procedures  specified  in 
Chapter  5.  Each  HOS  group  (Aviation,  Combat  support.  Combat 
Service  Support,  and  Ground  combat)  was  unigue  in  its  set  of 
functions.  However,  the  methcdology  was  similar  in  manipu¬ 
lating  the  data. 

a.  Function  BEFORE 

The  function  ”BEF0fiE**  performed  the  preliminary 
work  toL  the  parameters  tor  the  James-Stein,  Minimax,  and 
Maximum  Likelihood  estimation  processes.  See  Figure  0.2). 

b.  Function  GETGC81 

The  function  “GEIGCai”  derives  the  Normal  values 
for  the  James-Stein,  Hinimax,  and  Maximum  Likelihood  estima¬ 
tors.  See  Figure  D. 21. 
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• 

■"  —  '  —  - -  - - - -  ■  ■  - - - — 

7  BEF0RE:I5;I8 

fl  DOES  THE  PRELIMINARY  WORK  FOR  THE  PARAMETERS 

'2' 

A  FOR  THE  JAMES  STEIN,  MINMAX  AND  MAXLIKEHOOD 

■3 

A  FOR  THE  NORMALS  OF  THE  GROUND  COMBAT  MOS. 

■4 

I5^YG5  BINPREP  NG5 

'5' 

I8-YG8  BINPREP  NG8 

’6 

JG5-DLT  JAMES  15 

■7 

JG8-DLC  JAMES  18 

'8 

PG5«-JG5  BINCONV  NG5 

■9 

PG8^JG8  BINCONV  NG8 

’10' 

PM5*-YG5  MINMAX  NG5 

'll' 

PM8-YG8  MINMAX  NG8 

[12] 

IM5-(+/NG5)  ARCSIN  PM5 

‘13' 

IM8^(+/NG8)  ARCSIN  PM3 

'14' 

IA5-(+/YG5)  BINPREP(+/NG5) 

’15 

IA8^(+AG8)  BINPREP(+/NG8) 

'16' 

PA5-IA5  BINCONV  NC5 

[17‘ 

PA8<-IA3  BINCONV  NG8 

_  _  j 

Figure  i).20  API  Function  BEFORE. 


c.  Function  AFTER 


The  function  ••AFTER”  is  an  example  of  the 
conscliaation  of  th«  execution  of  multiple  years^  validation 
functions.  See  Figure  0.22. 


d.  Function  GETFOn 


The  function  •'GETFOfl”  derives  the  figure  of 
merit  (i.e.  mean  squared  plus  variance)  of  a  given  array  of 
values.  See  Figure  0.2J. 


Function  TEST 


The  function  "TEST''  was  utilized  in  the  addi¬ 
tional  comparison  suttests  mentioned  in  Chapter  5.  The 
function  conducted  the  two  comparison  tests  for  the 
projection  capability  cf  the  James-Stein  estimation  process. 
See  Figure  0.24. 


V  GETGC8 1  ;CS8 1  ;Y;N5;N8;Y5;Y8;IG5;IG8 
’1]  «  GETS  THE  NORMALS  FOR  THE  JAMES-STEIN,  MINMAX, 
’2]  «  AND  MAXLIKE  FOR  ThE  GROUND  COMBAT  MOS  FOR  THE 
'3]  «  PROJECTED  YEAR  1981.  ALSO  COMPUTES  THE  CHI 
’4]  ft  SQUARE  VALUES  FOR  VALIDATION  PURPOSES. 

’5]  l*-0 

'6]  N5**N81G5 

7]  N8-N81G8 

;8]  Y5«-Y81G5 

9]  Y8«-Y81G8 

'lO]  IG5^(Y5  BINPREP  N5) 

’ll]  IG8^(Y8  BINPREP  NS) 

’12]  NORaiG5-lG5-JG5 

'l3]  N0R81G8«-IG8-JG8 

’14]  DD*-DLT 

’15]  P<-PG5 

’16]  CS81G^(2  6)pO 

’17]  CS81G[1;  1  2]<^Y5  CHI  N5 

;i8]  DD-DLC 

19]  P^PG8 

;20]  CS81G[2:  1  2]‘-Y8  CHI  N8 
’21]  NMM81G5-IG5-IM5 
’22]  NMUaiGS^-lGa-IMS 

23]  DD-DLT 

24]  P^PMS 

’25]  CS81G[1;  3  4]-Y5  CHI  N5 
'26]  DD^DLC 
'27]  P-PM8 

’28]  CS81G[2;  3  4]-Y8  CHI  N8 
;29]  NASIGS-^IGS-IAS 
30]  NA81G8'^IG8-IA8 
'31]  DDi-DLT 
72]  P^PAS 

'33]  CS81G[1;  5  6]*-Y5  CHI  N5 
■34]  OD^DLC 
'35]  P»-PA8 

;36]  CS8IG[2;  5  6]^Y8  CHI  N8 

■37]  '  JAMES  OF  MINMAX  DF  MAXLIKE 

'38]  0-CS81G 


Figure  0.21  iPL  FuDctioA  GETGC81. 
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7  AFTER 

a  RUNS  THE  81  TO  83  PROJECTION  FUNCTIONS 
■GROUND  COMBAT  1981' 

GETGC81 

■GROUND  COMBAT  1982' 

GETGC82 

■GROUND  COMBAT  1983* 

GETGC83 


Figure  D.22  API  Function  AFTER. 


7  Z-GETFOM  X;M;N;V 

;i]  «  GETS  THE  FIGURE  OF  MERIT  (FOM)  OF  X 

\2]  fl  FOM  IS  THE  SQUARED  MEAN  PLUS  THE  VARIANCE 

;3]  M^(4/.X)-rN^(p.X) 

>]  V-(+/.((X-M)*2))4N-1 

5]  Z-V+(M*2) 

e]  ■MEAN  =  •.(»M) 

.[7]  'VARIANCE  =  '.(W) 

[3]  'FOM  =  '.(•Z) 


7  Z< 
[1] 
[2 
[3] 


[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 


19 


[20' 

[2i; 

[22; 


-Y  TEST  N;A;B;C;DD;IS:J;P;MLE;AC;SSEM;SEM;SSEP;SEP;DIFP;DIFM;NUM 
«  Y  IS  LEAVERS.  SHAPE  (2  YEARS.  31  LOS.  MOS'S) 

0  D  IS  THE  GLOBAL  INCIDENT  MATRIX  FOR  THE  DESIRED  GRADE 
•»  N  IS  CENTRAL  INVENTORY.  SHAPE  (2  YEARS,  31  LOS  MOS'S) 

IS-Y  8INPREP  N 

J^-D  JAMES  IS 

P-(J  BINCONV  N)xD 

MLE‘-((Y[l;;])^(N[l:;]+(NCl;;]=0)))xD 

AC-((+A)-^((+/N)+((+/N)=0)))xD 

SSEM*- +/.SEM  -  (AC  -  MLE)  .2 

SSEP-+/.SEP^(AC“P)*2 

DIFP-+/.(SEP<SEM) 

DIFM-+/.(SEM<SEP) 

NUM«-+/.D 

■1  THIS  IS  TEST  1.  ALL  THE  CELLS  IN  THE  FEASIBLE  REGION* 

•  ARE  ESTIMATED  BY  THE  JAMES  STEIN  PROCESS* 

■  SEP  =  THE  SQUARED  DIFFERENCE  OF  THE  JS  PROJECTED  AND  ACTUAL* 

*  SEM  =  THE  SQUARED  DIFFERENCE  OF  THE  MLE  AND  ACTUAL* 

•  SSE(MLE)  =  THE  SUM  OF  THE  SQUARED  DIFFERENCES  OF  SEM* 

*  SSE(P)  =  THE  SUM  OF  THE  SQUARED  DIFFERENCES  OF  SEP* 

I  « 

*FEASIBLE  CELLS  SEP<S£M  SEM<S£P  SSE(MLE)  SSE(P)* 

A*-‘  ' 


[23] 

[24] 

[25] 

[26] 
[27] 

28] 

’29] 

■30 

■31 

■32 

■33 

■34' 

'35' 

'36 

'37' 

38 

■39 

[40 

[4i: 

[42] 

[43] 
[44‘ 
[45^ 
[46' 


B^* 

C*-' 

DD«-'  ' 

B.(»NUM).A.(*DIFP),B.(*DIFM),C.(»SSEM),DD.(*SSEP) 

(  t 

SSEM^+/.SEM'-(AC~(MLEx(AC?^0)))*2 

SSEP«-+/.SEP»-(AC-(Px(aC^0)))*2 

DIFP4-+/,(SEP<SEM) 

DIFM-+/,(SEM<SEP) 

NUM-+/,(AC?^0) 

■  THIS  IS  TEST  2.  ONLY  THE  CELLS  IN  THE  ACTUAL  REGION  >0* 

■  ARE  REPRESENTED  IN  THE  COMPARISON  TEST  * 

■  SEP2  =  THE  SQUARED  DIFFERENCE  OF  THE  JS  PROJECTED  AND  ACTUAL* 

■  SEM2  =  THE  SQUARED  DIFFERENCE  OF  THE  MLE  AND  ACTUAL* 

*  SSE2(MLE)  =  THE  SUM  OF  THE  SQUARED  DIFFERENCES  OF  SEM2' 

*  SSE2(P)  =  THE  SUM  OF  THE  SQUARED  DIFFERENCES  OF  SEP2* 

(  I 

*FEASIBLE  CELLS  SEP2<SEM2  SEM2<SEP2  SSE2(MLE)  SSE2(P)* 

A*-'  ' 

B«-* 

C«-* 

DD-* 

B.(»NUM),A,(*DIFP),B,(»DIFM).C.(»SSEM).DD.(»SSEP) 


Piqure  D-24  IPI  Punctioa  TEST. 
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